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The  following  is  a brief  summary  of  the  work  performed  In 
this  contract  and  also  Indicates  certain  Important  open  problems 


More  details  can  be  found  in  the  attached  papers 


I.  SUMTIARY  OP  RESEARCH 


Our  work  was  motivated  by  the  potential  of  greatly 


improved  optical  communication  that  would  result  from  the  use 
of  a new  class  of  radiation  quantum  states.  These  states 
have,  been  named  two-photon  coherent  states  (TCS)^^\  and  a 
thorough  study  of  their  properties  has  been  completed^^^”^^ ^ . 

In  particular,  their  possible  generation  by  degenerate 
parametric  processes  and  two-photon  lasers  has  been  investigated 
Realization  of  the  quantum  measurement  that  leads  to  the 
predicted  quantum  noise  reduction  has  also  been  found  ' , 

The  point  process  counting  statistics  of  TCS’s  and  their 
performance  in  conventional  optical  reception  were  worked  out'-''' 
The  main  topic  of  this  contract  is  on  the  communication 
theoretic  aspects  of  TCS^ systems,  ^e-havm  shown  the  optimality 
of  TCS,  as  transmitter  quantum  states,  in  a number  of  linear 


modulation  and  binary  detection  problems 
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a theory  of  quantum  field  propagation 


which  Is  essential  for  a correct  analysis  of  TCS  behavior 
In  optical  transmission  systems.  In  particular,  we-ndw—" 
5iindersi;and  the  difference  between  near-fleld  and  far-fleld 

iS  c ?'  *» 

TCS  comr.iunicatlon.  While  one  can  realize  the  full  TCS  per- 


formance potential  in  near-field  communication,  the  far-fleld 
case  Is  uncertain  at  present. 

II.  SUGGESTIONS 

The  most  Important  open  problems  relating  to  TCS  communi- 
cation can  be  identified,  at  this  stage,  to  be  the  following: 

1.  Investigation  Into  the  posslbllty  of  employing 
TCS -state -generating  local  oscillators  for  greatly 
Improved  far-fleld  optical  communication. 

2,  To  design  and  subsequently  perform  an  experiment 
which  demonstrates  that  the  desired  quantiun  states, 

TCS’s,  can  be  produced. 

3.  To  ascertain  the  engineering  requirements  and 
sensitivity  tolerance  In  a practical  TCS  communi- 
cation system,  for  both  near  and  far  field  transmission. 

4,  To  develop  the  general  quantum  communication 
theory  of  optimum  state  selection  in  different 
communication  situations. 

III.  PUBLICATIONS  UNDER  THIS  CONTRACT 

1.  H.P,  Yuen,  "States  that  give  the  maximum  slgnal-to-nolse 
ration  for  a fixed  energy".  Physics  Letters,  vol,  56A,  pp.  105- 
106,  1976. 

2.  H.P.  Yuen,  "Tv/o-photon  coherent  states  of  the  radiation 
field".  Physical  Review  A,  vol.  13,  PP.  2226-2243,  1976. 
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3.  H.P.  Yuen,  ”Gausslan  states  and  parametric  processes,”  to 
be  submitted  to  Physical  Review. 

4.  H.P.  Yuen  and  J.H.  Shapiro,  "Optical  communication  with 
two-photon  coherent  states  - Part  I:  quantum  state  propaga- 
tion and  quantum  noise  reduction”,  to  be  published  In  IEEE 
Transactions  on  Information  Theory. 

5.  H.P.  Yuen  and  J.H.  Shapiro,  "Optical  communication  with 
two-photon  coherent  states  - Part  III:  quantum  measurements 
realizable  with  photoemlsslve  detectors",  to  be  submitted 

to  IEEE  Transactions  on  Information  Theory. 
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STATES  THAT  GIVE  THE  MAXIMUM  SIGNAL-TO-QUANTUM  NOISE  RATIO 

FOR  A FIXED  ENERGY' 
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Under  a radution  power  constraint,  the  maximum,  signal-to-quantum  noise  ratio  obtainable  for  any  state  of  a radi- 
ation field  is  found.  This  maximum  value  is  achieved  by  the  tw-o-photon  coherent  states  introduced  previously  to  de- 
scribe two-photon  lacers. 


j 

.1 


Let  the  density  operator  p be  an  arbitrary  quantum 
state  of  a radiation  mode  of  frequency  to  with  photon 
annihilation  operator  a.  We  write 

a=a^+ia2',  aj  =a|, 

for  selfadjoint  aj , 02-  The  signal -to-quantum  noise  ra- 
tio defined  by 

(S/AOp  s (tr  pay)'^ltr  p(Aaj )2  (2) 

with  Aaj  =aj  - tr  paj  is  a useful  measure  of  the  in- 
trinsic quantum  fluctuation  in  many  optical  experi- 
ments. In  particular,  if  the  operator  aj  is  measured  at 
a receiver  (optical  homodyne)  with  the  transmitted 
signal  in  state  p,  (S/N)^  provides  a measure  on  how 
well  the  signal  level  trpaj  can  be  estimated  from  the 
measurement.  We  first  show  that  under  the  energy  (or 
power  per  unit  frequency)  constraint 

ftco  tr  pa^a  < ftw  (3) 

for  a given  N^,  the  optimum  (S/N)p  achievable  by  any 
p is 

(S/N)^=4N,(,N,+  l).  (4) 

Note  that  under  (3),  trp(Aaj)^  cannot  be  made  arbi- 
trarily small  because  that  requires  an  arbitrarily  large 
as  apparent  in  the  following. 

Eq.  (4)  can  be  derived  by  using  the  uncertainty 
principle 

trp(Aa|)^  trp(Aa2)^  > 1/16.  (5) 

It  is  easy  to  show  that  all  the  available  energy  in  (3) 

* Work  supported  by  the  National  Aeronautics  and  Space 
Administration  under  Grant  NGL  22-009-01 3 and  by  the 
Office  of  Naval  Research  under  Contract  NR  042-342. 


should  be  used  and  that  no  energy  should  be  spent  in 
tr  pa2  so  that  (3)  becomes 

tr  p(Ajj  )2  ■¥  trp(Aa2)^  + (tr  paj  )2  = .iV^  + ^ . (6) 

From  (5)  it  is  clear  that  the  minimum  U p(Aa 2)^  re- 
quired to  satisfy  the  uncertainty  principle  for  a given 
tr  p(Aai  )^  is  [ 1 6 tr  pfAa,  )^  ] ' * . Thus  we  can  write 
(2)  in  terms  of  tr  p(Aa| )- alone  for  the  purpose  of  op- 
timization 

(S/7V)p  = {/V,+i-trp(Aai)2 

(7) 

— [16  trp(Aa|)2]  -t}/trp(Aai)2. 

The  maximum  value  of  (7)  is  readily  found  to  be  (4) 
with  corresponding 

tr  pai  = [N,iN,  + 1 )I(2A\  -i- 1 )]  • /2  (8) 

trp(Aa,)2=  (4(27Vs-H)]-l.  (9) 

This  optimum  value  is  achieved  by  the  two-photon 
coherent  states  |/3>g  previously  referred  as  “generalized 
coherent  states”  [1  ] . These  l/3>g  are  eigenstates  of  b 
= pa  + Mi'l  with  eigenvalues and  lp|2  — |i;|2  = l. 
When  p = 0,  |^>g  become  the  ordinary  coherent  states 
[2] . With 

0=  [/Vs(/V5  + 1)]1/2  (10) 

H = (N,+  \)I(2N,  + l)»/2,  u = NJ(.2N,  + l)>/2  (11) 

one  finds  from  ref.  [1]  that  g(/3|ail/3>  and 
g(/?|(Aa,)2l/3>g  are  indeed  given  by  (8)-(9). 

The  state  l^g  with  parameters  (lO)-(1 1)  is  an  or- 
dinary minimum  uncertainty  state.  In  general,  an  ar- 
bitrary two-photon  coherent  state  minimizes  the 
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ratio  (tr  pa\  )2/tr  p(A<i|  )^  for 

a\  =aj  cos^  + flj  sin0 

(12) 

a.  - . - v'p) 

0 = arc  tan — 

(13) 

2|p||p| 

The  maximum  value  (4)  obtained  from  1^)^  can  be 
compared  to  the  value 

(S/W)|„*,=  4rV,  (14) 

obtained  with  a coherent  state  p = laXal . By  spending 
a portion  of  available  energy  as  quantum  noise  in  02  to 
decrease  the  quantum  noise  in  aj , the  resulting  signal- 
to*quantum  noise.ratio  is  increased  from  a coherent 
state  as  if  JV^  is  effectively  increased  to  + JVy 
Among  other  things,  tliis  improvement  in  (S/yV)^  leads 
to  a higher  information  capacity  for  optical  communi- 
cations with  states  l|3>g  [3] . 

The  two-photon  coherent  states  1/3)^  can  be  pro- 
duced from  ideal  two-photon  lasers  [1 , 4]  or  degener- 


ate parametric  amplifiers  (3-6] . Not  only  do  they 
have  absolutely  small  <^af)  discussed  previously  [1,4], 
they  also  lead  to  higher  (5/A')  for  a fixed  radiation 
power  in  the  state  as  shown  above.  If  these  |/3>g  can 
be  experimentally  generated  from  two-photon  stimu- 
lated processes,  they  may  lead  to  significant  improve- 
nent  in  optical  system  performance.  They  may  also 
find  applications  in  physical  experiments  where  a (5/A') 
beyond  the  one-photon  laser  limit  (14)  is  sought  for 
low  noise  sensitivity 
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Two-photon  coherent  states  of  the  radiation  field* 
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The  concept  of  a iwo-pholon  coherent  sute  is  introduced  for  applications  in  quantum  optics.  It  is  a simple 
generalization  of  the  well-known  minimum-uncertainty  wave  packets.  The  detailed  properties  of  two-photon 
coherent  states  are  developed  and  distinguished  from  ordinary  coherent  states.  These  two-photon  coherent 
states  are  mathematically  generated  from  coherent  states  through  unitary  operators  associated  with  quadratic 
Hamiltonians.  Physically  they  are  the  radiation  states  ol'  ideal  two-photon  lasers  operating  far  above  threshold, 
according  to  the  self-consistent-field  approximation.  The  mean-square  quantum  noise  behavior  of  these  states, 
which  i'  oasically  the  same  as  those  of  minimum-uncertainty  states,  leads  to  applications  not  obtainable  from 
coherent  states  or  one-photon  lasers.  The  essential  behavior  of  two-photon  coherent  states  is  unchanged  by 
.nutl  losses  in  the  system.  The  counting  rates  or  distributions  these  sutes  generate  in  photocount  expenmentt 
also  reveal  their  difference  from  coherent  stales. 


I.  INTRODUCTION 

In  a recent  study*  of  quantum  communication  the- 
ory it  was  found  that  if  optical  radiations  can  be 
generated  in  a minimum-uncertainty  state**’  rather 
than  a coherent  state, ^ significant  performance 
improvement  may  result  in  an  optical  communica- 
tion system.  For  a fixed  radiation  mode  of  fre- 
quency it)  with  photon  annihilation  operator  a,  the 
coherent  states  | a>  (o  | a)  = a | a))  have  relatively 
large  quantum  fluctuations  when  ui/2ir  ^ 10*’  Hz. 
Writing  a for  self-adjoint  a,,a2,  a coherent 

sUte  |a)  gives  <Aa?>  =(AaJ>  = i.  A noise  energy’ 
is  then  obtained  if  either  a,  or  a,  is  measured 
in  optical  homodyne  detection.  This  quantum  noise 
is  frequently  dominant  over  other  noise  sources  in 
optical  communications;  for  example,  an  equiva- 
lent noise  temperature  of  a‘3400°K  is  obtained  for 
iJlw  at  the  YAIG  laser  frequency.  However,  states 
with  <AaJ>«  j and  correspondingly  larger  <Aa’> 
are  permitted  by  the  uncertainty  principle 

such  as  the  minimum  uncertainty 
states  with  <Aa’>  <Aaj>  = ji.  If  one  can  generate 
such  small  <Aa’)  states  by  an  explicit  physical 
process,  it  should  not  be  surprising  that  they  may 
be  profitably  used  in  many  applications. 

Minimum-uncertainty  states  cannot  be  obtained, 
however,  from  available  optical  sources,  all  of 
which  generate  coherent  states  and  their  random 
or  classical  superpositions.  An  investigation  of 
the  possible  ways  to  generate  minimum-uncer- 
tainty states  leads  to  a broad  class  of  radiation 
states  that  we  call  "two-photon  coherent  states” 
(TCS’s).  They  include  the  minimum-uncertainty 
states  of  which  the  coherent  state  is  a special 
case.  A detailed  study  of  these  TCS’s  is  reported 
in  this  article.  A few  of  cmr  results  have  been 
briefly  stated  in  Ref.  6. 


In  general,  a two-photon  coherent  state  differs 
from  a coherent  state  in  several  ways;  they  are 
generated  by  different  photon  processes,  they  have 
different  quantum  statistical  properties,  and  they 
have  different  coherence  properties.  Basically, 
coherent  states  are  generated  from  ideal  one- 
photon  stimulated  processes,  whereas  TCS’s  are 
obtained  from  ideal  stimulated  two-photon  process- 
es for  two  photons  of  the  same  mode.  The  usual 
multiphoton  parametric  processes’  are  one-photon 
processes  in  the  present  sense,  as  they  involve  a 
single  photon  from  each  separate  mode.  The 
quantum  noise  properties  of  TCS’s  are  basically 
the  same  as  those  of  minimum-uncertainty  states. 

Minimum-uncertainty  states  and  coherent  states 
can  be  easily  confused.  The  above  illustration  in 
terms  of  (AaJ)  and  is  a good  indication  of 

their  difference,  if  we  note  that  they  describe  the 
fluctuations  of  the  field  variables  directly  mea- 
sured in  a homodyne  experiment.  Further  careful 
discussions  will  be  found  in  Sec.  Ill  E and  else- 
where in  the  paper. 

Minimum-uncertainty  states  were  studied  at  the 
very  early  days  of  quantum  mechanics;  the  Ken- 
nard  packets'  are  particularly  familiar  in  the  con- 
text of  electrons.  It  is  therefore  no  surprise  that 
much  previous  work  has  touched  on  various  as- 
pects relating  to  the  present  two-photon  coherent 
states,  including  much  work  on  quantum  field  the- 
ory and  interacting  Bose  fluids.  It  is  impossible 
for  the  present  author  to  give  an  accurate  histori- 
cal survey  of  the  literature  pertaining  to  minimum- 
uncertainty  states  and  their  generation  by  quadrat- 
ic Hamiltonians.  The  particularly  relevant  papers 
now  known  to  the  author  include  Refs.  9-13  in 
quantum  optics  and  Refs.  14-17  in  quantum  field 
theory.  Nevertheless,  it  appears  that  TCS’s  have 
not  been  systematically  studied  before,  at  least  in 
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the  context  of  quantum  optics.  These  and  other  ad- 
ditional papers  will  be  mentioned  in  relation  to 
some  of  our  results  whenever  appropriate. 

This  work  is  not  at  all  a study  of  possible  math- 
ematical generalizations  of  coherent  states.  In 
fact,  the  TCS’s  can  be  considered  as  “coherent 
states”  in  a broad  sense,  as  discussed  later.  How- 
ever, they  should  be  distinguished  from  the  coher- 
ent states  of  Glauber,^  because  in  quantum  optics 
the  term  coherent  state  is  usually  understood  in 
the  Glauber  sense  which  carries  the  implication 
of  equal  fluctuations  in  the  quadratures  a,  and  a,. 

In  Sec.  H we  establish  notations  and  briefly  re- 
capitulate the  basic  results  in  quantum  optics  used 
in  this  paper.  For  a single  degree  of  freedom,  the. 
definition  and  properties  of  two-photon  coherent 
states  are  developed  in  Sec.  Ill,  together  with  a 
discussion  of  the  contrast  between  TCS’s  and  ordi- 
nary coherent  states.  The  generation  of  TCS’s 
from  quadratic  Hamiltonians  is  presented  in  Sec. 
rV,  with  illustrative  examples.  The  general  struc- 
ture and  properties  of  the  corresponding  unitary 
transformation  are  fully  determined.  Physical 
interpretation  of  quadratic  Hamiltonians  is  devel- 
oped in  Sec.  V in  terms  of  stimulated  two-phcton 
processes,  which  include  two-photon  lasers  and 
degenerate  parametric  amplifiers.  We  argue  that 
TCS’s  describe  ideal  two-photon  lasing  states  in 
analogy  to  the  coherent- state  description  of  one- 
photon  lasers.  The  approximation  which  leads  to 
this  conclusion  is  formalized  in  the  Appendix.  We 
also  show  that  the  mean- square  noise  behavior  of 
TCS’s  remains  essentially  the  same  when  the  ef- 
fect of  small  damping  and  additive  noise  are  in- 
cluded. In  Sec.  VI  we  briefly  discuss  the  coher- 
ence properties  and  counting  statistics  of  two-pho- 
ton coherent  states.  The  usefulness  of  TCS’s  in 
certain  applications  is  brought  out  in  Sec.  Vn.  The 
detailed  theory  of  multimode  TCS’s  will  be  given 
in  a separate  treatment. 

a COHERENT  STATES  AND  OPERATOR 
REPRESENTATION 

For  a single  mode  of  the  field  with  photon  anni- 
hilation operator  a,  a coherent  state*'' **  |a)  is  an 
eigenstate  of  a, 

a|a)«a|a).  (2.1) 

It  can  be  obtained  from  the  vacuum  state  |0)  via 
the  unitary  displacement  operator  D(o)»e“*^*“  *, 

|a)-D(a)|0).  (2.2) 

From  (2.2), 

<n|a)  = o"(«!)-'/V'*'».  (2.3) 

If  the  mode  has  a harmonic  time  variation  with 


frequency  u),  we  have  the  usual  canonical  variables 
q^(n/lwY>\a'  *a),  -a).  (2.4) 


In  general,  we  can  introduce  two  dimensionless 
self-adjoi.nt  operators  (a,,  a,). 


(2.5) 

} 

We  also  write 

a*a,+faj;  a„  a,  real. 

(2.6) 

and  refer  to  (a„aj)  or  (a„  a^)  as  the  quadrature 
components  (quantum  or  classical)  of  the  radiation 
mode.  These  quadrature  components  are  the  nat- 
ural variables  in  describing  a TCS,  even  when  the 
mode  has  a single  frequency.  Defining  the  eigen- 
states |a,),ofa„  {a,),  ofa,, 

a,|a,),  = a,|a,>„  a,|a2)j  = a2|a2>2.  (2-7) 

we  have 

i(aja')  = (2/ir)‘''exp[-(a,-a')*  + ia'(a'-  a'*)]. 

(2.8) 

From  (2.6)  it  follows  easily  that 

(a)  = a,  <JV)=|a|',  (2.9) 

(A«J>  = <Aa5)  = i,  (2.10) 

where  (AAf')»((Af-  (Af ))')  for  an  operator  M. 

The  angular  bracket  (M)  represents  the  expected 
value  trpM  of  M with  respect  to  a quantum  state 
p which  may  be  either  pure  or  mixed.  The  uncer- 
tainty principle  on  (a„a,)  is 

(2.11) 

Thus  I a)  is  a minimum-uncertainty  state.  From 
the  arithmetic- geometric  mean  inequality,  (2.11) 
implies 

(2.12) 

so  that  I a)  can  be  characterized  as  the  states  with 
both  minimum-uncertainty  product  and  minimum- 
uncertainty  sum. 

The  differential  operator  representation  of 
M(a*,a)  in  the  coherent- state  representation  is”’** 

<alM(«»,a)|)(i)»Af(a*,|  + ^)<a|<(.),  (2.13) 

from  which  one  can  show 

<a|Af(a’,fl)|8)*(al8)Af<’"(a*,8),  (2.14) 

where 

Af‘’”(a»,a)»<a|Af|a)  (2.15) 

is  the  normal-order  form’  of  M.  Note  that  it  is 
important  to  keep  the  separate  dependence  of  .If*"’ 
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on  a’  and  a in  (2.14)— A/*"’  is  a function  of  two 
comple.x  variables  (o’,  a).  Equation  (".14)  explic- 
itly shows  the  fact'*  that  a)  un.quely  de- 
termines A/.  In  terms  of  a;/,/j,),  the 

Schrodinger  equation  iHBU  '91  ^H(a' ,a)V  becomes'* 


iK 


at/'"* 

~a7~ 


=/f 


l/'"*  . 


(2.16) 


These  results  are  at  least  valid  for  bounded  oper- 
ators or  unbounded  operators  which  are  polynomi- 
als in  a and  a*. 

In  a number  of  mathematical  senses,*'***  any 
density  operator  p can  be  expressed  in  the  diago- 
nal  P representation,*^ 


p= J P(a",  o)|ot)(a|d*a,  d*a*d(Rea)d(Ima) . 

(2.17) 

Xj»(’)*f’1)  b®  ^be  normally  ordered  character- 
istic function 

(2.18) 

then  P(a*,  a)  is  the  two-dimensional  Fourier 
transform  of  XAf(n*i»?)i“'*‘ 


P(a*,a)=  J'e’*“***''XA,(»)*,i))d*Tj/jr*.  (2.19) 

It  has  not  been  demonstrated  that  P(a",  a)  provides 
a convenient  tool  in  actual  calculation  when  it  is 
not  sufficiently  well  behaved. 

Being  the  trace  of  the  product  of  a bounded  oper- 
ator and  a trace-class  operator,  the  character- 
istic functions  X/*(»?*,»?).  Xa(’1*,>7).  and  Xw(V',V)  are 
always  well  defined, 

(2.20) 

X»(n*,>?)»tr(pc-"****'’).  (2.21) 

The  Fourier  transform  of  Xw(V*tV)  is  ihe  Wigner 
distribution  and  that  of  Xa(V*,V)  is  p‘"’(a*,  o),**’*' 
p<">(at*,a)»(a|p|a>, 


JV(a\a)  = j’xr(»l*.»l)e'*“-’’“VT,/)r*,  (2.22) 

Xa(V.V)- / P'"*(a*,  a)e~*-“"Va/ir.  (2.23) 

The  characteristic  functions  are  related  by 
Xa(v*,  v)  =c'"’'*xw(»i*,  n) =«■'"*’' *xr(n*,  v) , 

(2.24) 

which  also  implies  relations  between  P(oi*,a), 
p‘**(af",  a),  and  W{a*,a). 

The  function  p'"’(a*,  a)  is  always  well  behaved, 
in  contrast  to  P(a*,  a),  and  is  also  quite  conve- 
nient for  moment  calculations.  Furthermore,  it 
specifies  o uniquely  as  in  (2.14)  and  is  the  prob- 


ability density  function**  describing  the  outcome 
statistics  of  a simultaneous  quantum  measurement 
of  two  noncommuting  observables**'”  (or  a mea- 
surement of  the  non-self-adjoint  operator  a,  opti- 
cal heterodyning).*'  We  will  use  p‘"*(a",  a)  or  the 
characteristic  functions  to  represent  a density 
operator,  particularly  because  P(a*,  a)  for  a TCS 
is  too  singular.  In  terms  of  p'"’,  we  have 

p = J^e"*'**®"*p*"’(a*,a  )«*"*%"  **  d*a  d*7j  /x* , (2.25) 

trpAf  = j lxJ(»I*,>l)]*Xa('7*,»l)®'’’‘'<^Vir,  (2.26) 

if  X *{tl * F tj)  = if (Afe** *•«"•*)  exists. 

III.  TWMHOTON  COHERENT  STATES 
A.  Definition 

For  a fixed  radiation  mode  with  photon  annihila- 
tion operator  a,  let 

(3.1) 

for  a pair  of  c numbers  p,  i/  obeying 

IpI*-  |>'|*  = l.  (3.2) 

It  follows  from  (3.2)  that 

[6,6*]  = 1.  (3.3) 

Any  transformation  b{a*,a)  which  leaves  the  com- 
mutator Invariant  as  in  (3.3)  is  called  a canonical 
transformation.  The  change  of  variables  from 
(a,  a*)  to  (6,6*)  according  to (3.1) and  (3. 2)  is  there- 
fore a linear  canonical  transformation.  A theorem 
of  Von  Neumann’*" '*  asserts  that  every  canonical 
transformation  can  be  represented  as  a unitary 
transformation,  i.e., 

b{a\a)  = Ualf  (3.4) 

for  a unitary  operator  U.  The  structure  and  physi- 
cal realization  of  the  that  leads  to  the  linear 
canonical  transformation  (3.1)  are  discussed  in 
Secs.  IV  and  V; 

tun-vJ . (3.5) 

The  commutator  (3.3)  or  Eq.  (3.4)  provides  6 

with  properties  exactly  similar  to  those  of  a.  One 
obtains  the  following  in  a way  identical  to  the  usual 
derivation  for  S and  a:’  Let  A*,  be  the  “quasiphoton” 
number  operator, 

N,^b^b  = Ui_NU\.  (3.6) 

Then  has  discrete  positive  eigenvalues  n,  with 
ground  state  |0^, 

1V,|0,>.0,  (3.7) 

|w,>  = l/Jm>.  (3.8) 
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Similar  to  |w>,  the  states  |»))^)  can  be  expressed  in 
the  form 

|m,>  = (6»)"(.«l)-‘'"|0,>.  (3.9) 

They  are  complete  orthonormal.  Furthermore,  b 
acts  as  the  lowering  operator  for  |»»,).  The  quasi- 
particle  interpretation  of  b and  |m,)  is  discussed 
in  Sec.  IV  A. 

The  TCS’s  |S>,  are  defined  to  be  the  eigenstates 
of  6 with  eigenvalues  0, 

= (3.10) 

The  n,  V dependence  of  a TCS  has  been  suppressed 
in  the  notation  ||3),.  When  necessary,  we  will  use 
the  more  complete  notation  |/3;  p,  v).  With  i/  = 0, 

|/3),  become  the  ordinary  coherent  states.  In  our  ‘ 
later  references  to  a TCS,  it  should  be  implicitly 


understood  that  u*0. 

From  (3.4),  there  is  an  independent  of  S 
which  gives 

l0>,  = Ui.|0)  (3.11) 

for  each  state  |S).  Therefore  in  parallel  with  |s) 
|S>,  = D,(0)|O>„  D,(S)  *«**’■•’*,  (3.12) 

j‘|S)„(^|d*S/ir=/,  (3.13) 

,0|^'>,  = exp(S*^'-i|<3|*-5l^'l')-  (3-14) 


If  the  notation  |S;  p,  v)  is  not  used,  the  same  p,  v 
is  understood  to  apply  to  all  of  the  in  an  equa- 
tion. 

From  (3.13),  any  state  vector  \ip)  can  be  ex- 
panded in  terms  of  |S),.  Similarly,  diagonal  op- 
erator representation  can  be  developed, 

Af=  J M“>(S»,S)|^)„<^|d"^/ir,  (3.15) 

for  the  antinormal  order  form  .Vf‘**(S*,S)  of  M.  In 
the  |S),  representation,  we  have  [similar  to  (2.21)] 

,<fi\M{b\  &)  U)  = Af  (3*, I ^ ,<3 1#  . (3. 16) 

Other  properties  associated  with  | a)  carry  over 
identically  to  |3),. 

The  results  of  this  subsection  on  |3),  and  |m,) 
depends  only  on  the  relation  (3.4),  and  not  on  the 
more  specific  (3.5).  The  specific  characteristics 
of  (3.1)  are  developed  below.  Certain  properties  of 
TCS’s  have  been  discussed  by  Stoler’  in  the  con- 
text of  minimum-uncertainty  states. 

B.  Wave  functions 

The  general  coherent- state  wave  function 
of  a TCS  is  first  determined  as  follows:  From 
(3.1),  (3.10),  and  (2.13)  we  obtain 


The  solution  of  (3.17)  is  of  the  form 

<o|3>,  = exp[Q!*3/p-  2 |a|*-  i/ar*V2p4-/(3*,3)]. 

(3.18) 

The  (o*,  a)- independent  / can  be  determined  from 
(3.16),  (3.1),  (2.13),  and 

= (pa+^  + i'^),<3|a),  (3.19) 

so  that  the  3-dependent  part  of / in  (3.18)  is  -5  |3|* 
+ i/*3*/2p.  From  the  normalization 

J|<o|3>,|’</Wir  = l, 

<a|«a  is  thus  determined  up  to  a (a*,  a;3»,3)- 
independent  phase  in  the  form 

<«|0>a  = M‘‘'’exp[-.i|a|*-2|3|^-(v/2p)a** 

+(v/2p)3'+(l/M)a*0  + teo]  (3.20) 

for  a real  0,.  The  normalization  constant  is  cal- 
culated from  the  normalization  of  Gaussian  inte- 
grals.“ For  p = l,  »'  = 0,  (3.20)  reduces  to  the 
usual  form  of  (a|3)  for  0o  = ^>  ^ choice  of  phase 
we  will  adopt  for  (a  1 3),. 

While  we  have 

K|0>,  = 8’"(ml)-‘/V'»'’'S  (3.21) 

similar  to  (2.3),  the  function  (w|3),  takes  a more 
complicated  form  which  can  be  determined  as 
follows:  As  a function  of  the  complex  variables 
z and  t,  «***•*’  can  be  expanded  in  the  Taylor 
series 

kl<-.  (3.22) 

So 

where  /f„(z)  is  the  nth  Hermite  polynomial  with 
complex  argument  z.  Writing 

and  usir.g  (3.20)  and  (2.3),  we  find 

<n  |3>,  = (nl  p)-‘'*(i;/2p)"'X[^(2pv)-'*] 

xexp[-||3i*  + (i'V2M)3l.  (3.23) 

For  )i  = l,  u = Q,  the  asymptotic  forms  of  /f„(z)  for 
various”  argz  all  lead  to  (2.3),  as  they  must. 
Equation  (3.23)  shows  that  the  counting  statistics 
of  a TCS  is  far  from  the  Poisson  (2.3). 


iMHSB 


\ 


i 
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The  wave  function  can  be  readily  determined  from  (3.20)  and  (2.8).  Thus  we  find 


2\'/« 


(ti  _ »<)■*'*  exp 


(-* 


2(n-^)'  -M- 

From  (3.20)  we  have 

H,v\0a>  Mu.  Mo)  =(MoM*-  MoM*)"''’ 

xpxnf  Mal»  i|a|3,  ^ fl.fl  ■^oM-Mpt'  ...  WoM-MoM)*  A 

i\p\  i\Pa\  “ 2((XoM*- VoV»)^  2((ioP*  - v^V) 
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(3.24) 


(3.25) 


Note  the  Gaussian  form  of  (3.24)  and  (3.25).  Var- 
ious reduction  of  (3.24)  and  (3.25)  to  the  previous 
formulas  can  be  easily  observed.  The  relative 
phases  of  (3.20),  (3.21),  (3.24),  and  (3.25)  are 
mutually  consistent,  as  well  as  being  consistent 
with  the  above  wave  functions  for  the  coherent 
states. 

From 

53  <«  I I ^># 

for  the  quasiphoton  number  states  | w,)  we  obtain 
(«lw,>  = Qi(M,M:»i,n), 

///««,  w odd,  /{Odd,  (3.26a) 

= Qi(m,m;  w,n)+Q2(M,v;m,n), 

m even,  n even,  (3.26b) 

= Q2(m,«';  »»,«), 

ni>n,  ni  even,  n even,  (3.26c) 

= 0,  otherwise;  (3.26d) 

(_  2)‘"-">'*(n!)‘'’ 


Q,(M,  v;  »i,n)=p.' 


[(«-w)/2]!(mI)‘'"’ 
(3.26e) 


(nlm!)*''^ 

(n/2)!(»w/2)r 


(3.26f) 

Note  that  m and  rt  must  be  both  even  or  both  odd  in 
order  that  (n|»tP  be  nonvanishing.  Thus  only  an 
even  number  of  photons  would  be  counted  for  the 
state  |op=  |0),.  Other  wave  functions  of  Im,)  can 
be  obtained  accordingly.  These  wave  functions 
provide  the  statistics  of  various  measurements  on 
the  radiation  states  Ifl),  and  |/«,). 

C.  Quantum  fluctuations  and  characteristic  functions 
From  (3.1)  and  (3.2), 

a=p*ft-v6’,  (3.27) 

so  that  for  the  TCS  ]0),, 


(a>3,03la|6),=  M*6-v6*s6*6.  + i6j,  (3.28) 
<aV=  131"+ |«'|^  (3.29) 

(Aa^>  = - 4*v  = (Aa’»>»,  (3.30) 

(Aaf)  = j I M - mI^,  (3.31a) 


<l<7f>  = i|p  + v| 


(3.31o) 


(Aa,Aaj>5((a,-3i)(a2-  32)>  = ii(M*v-  +1) , 

(3.32a) 

(AnjArt()  = i((p*i'- I) . (3.32b) 

Equation  (3.32)  is,  of  course,  consistent  with  the 
commutator 

[a„aj]  = [6„6j]=i! . (3.33) 

Compared  to  the  case  of  a coherent  state  where 
(a|(A{i)"(Aa)|a)  = 0,  (3.34) 

a TCS  (/J),  has  added  quantum  noise  energy 

,O|(Aa)»(Aa)l0>,=  IH".  (3.35) 

For  self-adjoint  hj,  and  realSi.  Sz,  let 

6s6,  + i6j,  + (3.36) 

The  TCS  1 3),  can  be  characterized  as  the  states 
which  minimize  both  (Ah,)  (A6*)  » ^ and  (A6*) 

+ (A6*)  On  the  other  hand,  (3),  minimizes  the 
uncertainty  product  only  when 

p = 6v,  for  6 real.  (3.37) 

These  minimum  uncertainty  states’  include  the 
eigenstates  of  fli.Oj  in  the  limit  p.,  v-«,  such  that 
(3.2)  holds.  Specifically,  the  eigenstates  |a,), 
are  obtained  when  6-1  and  the  eigenstates  |aj>j 
are  obtained  when  6 — - 1.  As  with  the  minimum- 
uncertainty  states,  the  mean-square  quantum  fluc- 
tuations in  <1,  and  a,  be  exchanged  in  (3',  from 
(3.31).  On  the  other  hand,  they  are  fi-xed  at  } (or 
coherent  states. 

Thus  for  a mode  of  frequency  u>,  a significant 
noise  energy’  iffu)  is  obtained  if  either  of  the  qua- 
drature components  is  measured  in  a coherent 
state,  while  (Aap«  t with  corresponding  large  <Aa‘) 
can  be  obtained  for  measurement  of  a,  in  a TCS. 
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A mathematical  scaling  of  a will,  of  course,  lead 
to  i even  in  an  ordinary  coherent  state. 

For  example,  j is  possible  with  b given  by 

(3.1).  In  actual  applications  such  scaling  cannot 
affect  the  system  performance.  In  Sec.  VII B we 
will  show  how  the  noise  in  a coherent  signal  can 
be  reduced  if  |S),  is  available,  which  is  impossible 
by  scaling  on  |S).  The  fact  is  that  (Aa‘)  measures 
the  absolute  amount  of  noise  in  a quadrature  com- 
ponent and  a is  a more  fundamental  variable  than 


A a/  \ 

\ -2C,  1+2C,/ 


(3.47) 

(3.48) 


In  terms  of  the  variables  (a*,  a),  (3.38)  becomes 
p‘"’(a*,  a) » 1 p 1''  exp[-i(a,  - a,)’'A*‘(a,  - 5,)] 


cussion  and  clarification. 

From  (3.20)  one  readily  obtains 

through  the  transformation 

(3.49) 

l<a|^>,P=|p|-‘exp[-(l-2C.)(a,-3,)» 

-(U2C,)(a,-32)’ 

'3.50) 

+ 4C,(o.-3,)(a2-^2)l,  (3.38) 

where 

A..f 

<AnAa’>\ 

(Aa^*>  / 

(3.51) 

Ca_  i//2paC,+fC„  C„C,  real,  (3.39) 

and  ^ is  given  by  (3.28).  From  (3.2),  so 

that 


where  T denotes  the  transpose  of  a matrix. 

The  antinormally  ordered  characteristic  func- 
tlon  xi(»I*,»l)  for 


(3.40) 

Thus  (3.38)  is  indeed  a Gaussian  joint  probability 
density  function  in  (a,,  a,),  with  mean  and 

variance 

o?  aE[(a,  - 3,)’]  = id  + 2C,)/(1  - 4 Ic  I *) , (3.41a) 

(j*aE[(a,  - 3,)’]  = i(l  - 2Ci)/(l  - 4 |C  I’) . (3.41b) 

The  notation  E[R]  indicates  the  probabilistic  aver- 
age  of  a random  variable  R with  respect  to  its 
density  function.”  The  correlation  coefficient  is 

r s£[(at,  _ 3,)(aj-  = 2Cj(l  - 4C  ,)*''* . 

(3.42) 

This  density  function  is  illustrated  in  Fig.  1. 

Equation  (3.38)  is  conveniently  expressed  in  the 
matrix  form 

<*2)  ^ Ka 

* 1 p I*' exp{-i(a, - o,)’'A'‘(o, - a,)] , 


(3.43) 

where 

(3.44) 

5,  a £[£,]»  3,, 

(3.45) 

(3.46) 

FIG.  1.  Representation  of  a two-photon  coherent 
state  by  p^"*  (or  1,0 j)  from  (3.38),  with  of>ofi  a coherent 
state  ts  obtained  for  of  = 05 . 
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P,*|0>„<i3|  (3.52) 

is  then 

XiK'?!.  Vi)  = expUrj^a,  - , (3.53) 

(3.54) 

V-2T,,/ 


and  the  Weyl  characteristic  function  is,  from 
(2.32), 

y’wiVi,  Jfj)  = exot/^n,  - ■ (3.55) 

Since  A>0  (positive  definite), “ the  Wigner  distri- 
bution is  also  a Gaussian  density  function.  The 
antinormally  ordered  characteristic  function 
X*(»7,,»?a)  is 

Xjr(»Ji.  tJa)  =exp[inja,-  i5f(A  - Uju,]  ■ (3.56) 

The  antinormal-,  symmetrized-,  and  normal- 
order  moments  of  (a*,  a)  can  be  calculated  from 
(3.53),  (3.55),  and  (3.56),  respectively.  Note  that 
these  X)f>  Xw.  Xa>  p'"*  TCS’s  are  all 

Gaussia.i.  The  above  matrix  forms  are  convenient 
for  generalization  to  the  multimode  situations. 

The  minimum  noise  (A6;)  obtainable  in  a TCS 
for  a self-adjoint 

yj+y|=l,  (3.57) 

is  the  minimum  eigenvalue  of  A,,  with  the  mini- 
mizing value  of  (yi.yj)  given  by  the  corresponding 
eigenvector.”  The  result  can  be  expressed  in  the 
following  form:  Let 

(3.58) 

for  a real  ^ and  let 

o'aae**.  (3.59) 

Then  a\  minimizes  (Ab’)  when 

0atan**[»(p*i'-»'*p)/(2|p|  |i'i  * p*u*v*ti)] 

(3.60) 


and  the  minimum  fluctuation  is 

<Aa'»)=i(lp|-lvl)>.  (3.6la) 

The  conjugate  variable  is  the  one  which  max- 
imizes (Ab,),  with 

<A«'*)  = i(|MMvl)>.  (3.61b) 

From  (3.61),  the  TCS  |/3>,  can  be  seen  to  mini- 
mize the  uncertainty  product 

(Afl;»>(A«5»)si,  (3.62) 

for  a'  defined  by  (3.59)  and  (3.60).  Thus  |0),  are 
the  minimum-uncertainty-product  states  for  a' 
and  can  therefore  be  understood  as  generalized 
or  rotated  minimum-uncertainty  states.  This 


fact  is  important  in  some  applications  as  indicated 
in  Sec.  VII B,  because  the  original  (Ao^)  and  <Aa|) 
can  be  both  larger  than  } from  (3.31)  unless  p/v 
is  real. 

D.  Diagonal  operator  representation  for  I 

Let  0,  represent  a density  operator  which  can  be 
written  as  a random  superposition  of  coherent 
states,  i.e., 

o, =y  P.(a»,  a)  I a>  (Of  I d*a  , (3.63) 

for  a positive  true  probability  density  function 
P,(a*,a).  It  follows  that 

tro, (Aa?)  = i+ f P,(a„  aj)(a,  - a^fda^da^  s i , 

(3.64a) 

for 

“i=/ P,(c^^,a^)a^da^da^, 
and  similarly, 

trp, (Aa^)  » i . {3.64b) 

Equation  (3.64)  states  the  obvious  physical  fact 
that  the  introduction  of  classical  randomness 
merely  increases  the  noise  in  a,  and  a,. 

A TCS  with  (3.31)  therefore  cannot  be  a mixture 
of  cohere.nt  states.  In  fact,  the  density  operator 
p,  cannot  be  written  in  the  diagonal  form  (2.17) 
with  P{a*,  a)  belonging  to  the  class  of  tempered 
distribution.  This  can  be  seen  from  the  normally 
ordered  characteristic  function  Xj»(7ii  72)-  for  a 
TCS,  A - cannot  be  positive  definite.  Thus 
Xy  does  not  possess  a Fourier  transform.  If  we 
insist  on  having  a diagonal  P representation  for 
p,,  we  will  need  the  Fourier  transform  of  e”*  for 
a real  variable  x and  y>0.  It  does  not  appear  to 
be  useful  to  talk  about  the  Fourier  transform  of 
such  a quantity,  which  certainly  does  not  belong  to  r 
the  class  of  tempered  distribution.  Further  dis- 
cussion of  this  divergence  is  given  in  Ref.  27. 

One  may  still  consider  the  diagonal  P represen- 
tation of  p,  as  belonging  to  the  class  of  distribution 
or  alternatively  as  the  limit  of  a sequence 
of  well-behaved  density  functions,^^  Since  no  use- 
ful method  of  calculation  has  resulted  from  such 
viewpoints  and  since  o*’'’(a*,  a)  is  sufficient  for 
our  purpose,  we  will  say  that  the  P representation 
of  p,  does  not  exist  and  we  will  not  employ  it  in  the 
following. 

Note,  however,  that  a generalized  P,  represen- 
tation of  Of  does  exist  and  can  be  useful.  Thus 

0-f  Pf(fi\0)\^,.<fi\dW»  (3.65) 

may  be  used  to  represent  |3),,  with 


’fniiiiWi'liiU  r Maw 
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= (3-6®) 

Equation  (3.66)  also  describes  random  mixtures  of 
Is),  in  the  same  way  as  (3.63).  In  particular,  it 
remains  useful  with  the  introduction  of  an  additive 
noise.  However,  it  is  inconvenient  if  (p,  u)  are 
time  varying.  While  it  can  be  used  in  moment 
calculations  by  expressing  M{a^,a)  in  normal- 
order  form  in  terms  of  the  variables 

(6*,b),  it  is  no  longer  particularly  suited  to  be 
calculations  of  [a*,a)  normal-order  moments  en- 
countered in  the  description  of  most  experimental 
results.  Therefore  (3.65)  will  also  not  be  em- 
ployed in  this  paper. 

E.  Coherent  sutes  venus  two-photon  coherent  states 

As  mentioned  in  Sec.  I,  two-photon  coherent 
states  were  encountered  in  many  different  stud- 
ies, but  they  received  little  attention  in  quantum 
optics.  Obviously  this  is  due  largely  to  the  fact 
that  these  states  are  not  generated  from  available 
optical  sources,  but  other  factors  may  have  also 
contributed. 

It  is  well  known  that  the  mathematical  form  of  a 
minimum-uncertainty  state  at  a frequency  is 
identical  to  that  of  a coherent  state  at  a different 
frequency  w,.  Therefore  one  may  tend  to  view 
minimum-uncertainty  states  as  just  coherent 
states  of  a different  frequency,  particularly  since  ^ 
the  frequency  is  often  normalized  away.  This 
viewpoint  is  erroneous,  because  in  the  state  space 
3C(u),)®3C(w,)  a minimum-uncertainty  state  at  w, 
is  given  by  l0),®  |0),  while  a coherent  state  at 
is  (0)®|^).  Physically,  no  photocounts  would  be 
registered  for  a coherent  state  at  to,  in  a counting 
experiment  employing  detectors  responding  only 
to  fields  at  w,.  In  this  connection  one  may  also 
note  that  lor  a fixed  mode  the  frequency  is  also 
fixed  by  the  dispersion  relation.  Changing  the 
frequency  of  a fixed  mode  by  dielectric  modulation 
would  not  yield  a minimum -uncertainty  state  at  fre- 
quency Wj  from  a coherent  state  at  frequency  iOp  A 
coherent  state  at  a different  frequency  u),  would  be 
obtained  instead. 

Alternatively,  one  may  regard  TCS’s  and  ordi- 
nary coherent  states  as  (mathematically)  equiva- 
lent, since  they  are  related  through  a unitary 
transformation  or  some  kind  of  scaling.  The  ques- 
tion of  scaling  was  already  dealt  with  in  Sec.  IIIC 
concerning  the  noise  in  quadrature  components. 

It  is  true  that  two-photon  coherent  states  are  co- 
herent states  in  a broad  sense:  they  are  the  co- 
herent states  with  respect  to  b{v*0)  instead  of  a. 
However,  this  abstract  unitary  equivalence  means 
little  in  describing  actual  experimental  results. 
Thus  («,la)  =<M|a()  but  (njo),  has  little  resem- 
blance to  (»i  la).  In  general,  for  an  arbitrary  uni- 


tary transformation  U we  can  define  a set  of  states 

l>. 

|)  = tf|o).  (3.67) 

The  functions  (x  | ) cannot  be  obtained  by  simple 
scaling  on  (x  ( a)  and  has  to  be  computed  anew  via 
(x|Cf|a).  This  paper  develops  such  specific  re- 
sults for  the  Uj.  of  (3.5). 

With  respect  to  a,  whose  zero  eigenstate  is  the 
radiation  ground  state  occurring  in  nature,  the 
properties  of  ] o),  are  quite  different  from  | a)  in 
many  respects,  as  presented  in  this  section.  Since 
the  term  coherent  state  is  usually  understood  in 
the  original  sense  (2.1),  a new  terminology  is  re- 
quired for  finer  distinction  between  coherent  states 
in  the  broad  sense.  The  term  two-photon  coherent 
states”  is  adopted  because  of  the  close  connection 
between  these  states  and  quadratic  Hamiltonians 
or  two-photon  lasing  processes,  discussed  below. 

rv.  QUADRATIC  HAMILTONIAN 
AND  STRUCTURE  OF 

From  (3.11)  and  (3.20)  we  can  see  that 
is  an  exponentiation  of  a quadratic  form  in  (a*, a). 
Such  a Cf£"’(o*,a)  can  be  written  as  <ole“lo)  for  a 
self-adjoint  operator  Af  quadratic  in  (a^,a),  from  a 
theorem  of  McCoy Therefore  U j must  be  the 
unitary  operator  generated  by  a quadratic  Hamil- 
tonian Hq.  The  mathematical  solution  </g(i , i,)  of 
the  Schriidinger  equation 

ft  rr 

£/o«o,0  = l,  (4.1) 

is  investigated  in  some  detail  in  this  section,  which 
gives  all  of  the  Ui_. 

A.  Positive  definite  quadratic  Hamiltonian 

Consider  the  general  quadratic  Hamiltonian 

/f  0 = » UiO*a  +/ia  * ’ +/ * a +/3a’),  (4 .2) 

where  the  c numbers  /,  may  be  time  dependent. 
From  Hermiticity  /,  has  to  be  real.  For  a mode  of 
frequency  u,  we  have  /i  = «,  and  Kf^a*a  represents 
the  free  radiation  energy  of  that  mode.  The  /,  and 
/,  terms  then  represent  interaction  energies.  In 
particular,  the /,  terms  describe  a two-photon 
mechanism,  whereas  the  /,  terms  describe  the 
usual  one-photon  or  linear  driving  mechanism. 

If  a physical  system  with  Hamiltonic"  given  by 
(4.2)  is  stable,  Hq  should  not  have  a spectrum 
ranging  through  the  entire  continuum  from  to 
«o  in  order  for  a ground  state  to  be  defined,  i.e., 
the  spectrum  of  Hq  should  be  lower  bounded.  By  a 
shift  of  seal-*,  one  may  require  that  the  Hamilton- 
ian be  posi..ue  definite.”  If  this  condition  is  not 
satisfied,  (4.2)  cannot  be  properly  used  to  de- 
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scribe  a physical  system  of  interest  in  our  present 
context.  It  has  to  be  extended  to  Incluae  a more 
complete  description  of  other  relevant  .'actors  of 
Importance. 

Under  the  condition 


A>2|Al 

(4.3) 

the  Hamiltonian  (4.2)  can  be  changed,  similar  to 
the  Bogoliubov  transformation,’’  to  the  following 
form: 

Ro=/ob’b+/o(l4/,*-*'/al'+l‘'|-). 

(4.4) 

b*4a+vo^  + 4/J- v/j. 

(4.5) 

4*[2//o(A-/o)f'’/*e'*, 

(4.6a) 

«'»[(A-/o)//of 

(4.6b) 

(4.7) 

where  <>  is  an  arbitrary  phase.  The  canonical 
transformation  (4.5)  is  an  inhomogeneous  linear 
one,  with  (4.6)  obeying  (3.2).  The  Hamiltonian 
(4.4)  is  manifestly  positive  definite,  and  (4.3)  is 
indeed  the  necessary  and  sufficient  condition  for  a 
classical  Hq  of  the  form  (4.2)  to  be  positive  defin- 
ite. 

Similar  to  the  situation  of  an  interacting  Bose 
fluid,  h*  of  (4.5)  can  be  interpreted  as  the  quasi- 
boson creation  operator  with  a quasiparticle  spec- 
trum /q.  The  quasiparticle  ground  state  is  the  |o,) 
of  (3.7),  whereas  |»i,)  are  the  excited  states.  It 
seems  highly  unlikely,  however,  that  the  radiation 
states  I nr,)  may  be  generated  optically. 

The  following  observation  should  be  made;  A 
linear  canonical  transformation  from  a to  b can  al- 
ways produce  /,  terms  in  b from  a free  Hamilton- 
ian hf^a^a.  In  fact,  it  is  easily  checked  that  an  /, 
term  results  if  and  only  if  the  v term  mixing  b^  is 
present  in  the  transformation.  The  transformed 
Hamiltonian  contains  a constant  term  |ir|^,  so  that 
the  eigenstates  of  b have  this  added  excitation  ener- 
gy compared  to  the  true  ground  state  |0).  Such  a 
Hamiltonian  is  clearly  not  physically  equivalent  to 
(4.2).  In  this  sense  the  operator  a is  more  funda- 
mental than  b and  further  justifies  the  need  for 
viewing  |a),  as  a two-photon  coherent  state  with 
respect  to  a. 

Our  following  results  on  are  actually 

valid  independent  of  the  condition  (4.3).  Neverthe- 
less, one  should  note  the  significance  of  the  posi- 
tive definiteness  of  (4.2). 


B.  Explicit  solution  for  Vq  I/./o  ) 

The  solution  can  be  obtained  by  normal- 

ordering techniques,’  among  many  possible  meth- 
ods. With  (4,2),  the  solution  U{,">(a*,a;f,fo)  of 
(2.16)  is  of  the  form 


U<,">(a*,a;/,fo)'exp[A(/)+B(f)a’  + C{f)a‘’ 

+B(/)|al’+£(/)as-f(/)a*], 

(4.8) 

where 

^ =-  2(2/, T +/?£’+/,*£),  (4.9a) 

^=-f/?(D+l)",  (4.9b) 

^ = - 2(4/ ‘C*  + 2/,C  +/,),  (4.9c) 

^ = - i(Af*C  +/,)(£  + !),  (4  .W) 

^ =-  i(2/?F +/?)(£ + 1),  (4.9e) 

^ * [(4/?C  + A) F + 2/,*C  +A ] , (4.9f ) 


with  the  initial  condition 

A(/o)=B«o)=C(fo)=l>«o)=£«o)=f«o)=0-  (4-10) 
From  (4.8)  and  (2.14)  we  get 

<a|C/'g(<,/o)|/3)  = exp{-  i|a|’-  ^|i3|’+i4(f)+B(<)/3’ 

+ C(f)a*»  + [U(f)  + lla*^ 

+ £(f)0+F(f)a*}.  (4.11) 

Equation  (4.11)  can  be  written  in  the  form  of 
<ali3  + t(f);p(/)v(/)>,  apart  from  a phase  factor. 
Using  (3.20),  we  find  by  comparison  that 

Uo{t,0\p)‘ \l3*C{ty.n{t),u(t))e'»^‘\  (4.12) 

0(f)^  e(/)  + 1 i[/3*t(0  - ^£*(f)),  (4.13) 

for  a real- valued  (^,2;)  independent  function  6(f) 
and 

A(t)  = -i  Inp(f)  + v*(f)t’(f)/2p(f) 

-i!£(0i'  + ief),  (4.14a) 

B(f)  = v*(f)/2p(f),  (4.14b) 

C(f)*-i;(/)/2M(/),  B(/)  = l/#i(f)-l,  (4.14c) 

£(/)  = v*(/)»x(f)/£(f)  - £*(f),  £(f)  = £(f)/4(f). 

(4.14d) 

If  (4.14)  is  substituted  into  (4.9),  the  equations  of 
motion  for  4,  v,  andOare  obtained, 

^=iA4-2«/Jv,  ^=-rAv  + 2iAp,  (4.15) 

^ =f  (/,*•' -/,4), 
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Equations  (4.15)-(4.17)are  equivalent  to  (4.9).  The 
initial  condition  (4.10)  becomes 


M0)  = l,  i'(0)=0,  (4.18) 

£(0)  = 0,  (4.19) 

while 

0(0)  = 0 

is  obtained  from  (4.12).  The  general  solution  (4.8) 
or  (4.12)  can  be  given  in  quadrature  once  the  solu- 
tion for  (4.15)  is  given  in  quadrature,  since  the 
other  equations  can  then  be  immediately  integrated 
Defining 


•'‘o 


(4.21) 


(4.22) 


the  solution  of  the  first-order  coupled  differential 
equation  (4.15)  under  (4.18)  is  equivalent  to  the 
solution  of  the  following  equation  for  a new  vari- 
able e : 


^=/*e»-/.  e(0)=o. 


(4.23) 


es2Ce»‘''‘  ‘o>.  (4.24)  , 

Once  C is  obtained  through  (4.23)  and  (4.24),  all  of 
the  other  variables  A,  B , D,  E , and  F can  be 
given  in  quadratures  from  (4.19).  But  Eq.  (4.23), 
a case  of  the  well-studied  Riccati  equation,  does 
not  admit  a solution  in  quadrature  in  general,  al- 
though many  methods  and  results  on  its  solution 
are  available.'"''** 

The  existence  and  uniqueness  of  a solution  to 
(4.23)  imply  that  the  solution  of  (4.9)  and  (4.10)  al- 
so exists  and  is  unique.  Thus  (4.8)  is  indeed  the 
solution  we  seek.  With  tf*''’(a*,a)  given  by  (4.8), 
the  explicit  operator  expression  of  Ug(t,ta)  is 

U{t  t ) = • 

(4.25) 

Equation  (4.25)  can  also  be  written  as 
exp(y  js” + + Vafl  * + + >50) 
for  a set  of  c numbers  y,.'’ 

C.  Examples 

Equ.'ition  (4.23)  can  be  integrated  when  ///*  is 
independent  of  /,  i.e.,  when 

/a(0  = r(f)exp[i(|>,-  2iSr(f,lo)),  (4.25) 

(5j  independent  of  I,  for  a real- valued  function  r{f) 


which  can  take  on  both  positive  and  negative  values. 
Under  (4.26),  one  finds  readily  from  (4.23),  (4.24), 
(4.9d),  (4.14c),  and  (4.18)  that 


p(/)=e‘®'‘'‘o>cosh[2(R(/)l, 

!<(/)  = fe '•‘o’ sinh[2(R(/)], 

<R(f)»  r(y)dT, 


(4.20) 

...  -'<« 


(4.27a) 

(4J7b) 

(4.28) 

(4.29) 


Equation  (4.26)  is  especially  interesting  because 
of  the  resonant  factor  obtained  when  /, 

- (>).  This  factor  is  automatically  present  in  our 
later  (diysical  interpretation  of  (4.2).  The  noise 
behavior  in  this  case  is  of  great  interest.  Thus 
when  (with  also  /i  = w,  = 

= (4.30) 

r = const  independent  of  /,  we  have  from  (3.61)  and 
(4.27) 


(A«J»(<)>  = ie-*’-‘,  (An«(f))  = ie*’-‘,  (4.31) 

for  a'  =oe'"*.  In  this  case  the  noise  in  one  quadra- 
ture component,  apart  from  the  optical  oscillation, 
decreases  exponentially  while  the  conjugate  com- 
ponent has  exponentially  increasing  noise.  The 
quadrature  components  a^  and  a,  are  completely 
symmetrical  here.  We  can  choose  the  sign  of  r, 
or  equivalently  the  phase  of  r,  to  cause  either 
one  of  them  to  have  decreasing  noise.  For 
<l)j  = 0,  >'  = const,  we  again  have  (4.31)  with  a' 

= The  nature  of  phase  control  in  a'  from 

the  phase  of  should  be  clear  in  this  example. 

Equation  (4.15)  can  be  integrated  for  a constant 
/,  independent  of  time,  with 

p(f)  = cos(Xf)  + (fw/X)  sin(X/), 

i/(f)  = f(2/j/x)sin(x/), 

for  Xs(u>’-4l/jl*)*^*^0.  The  quantum  noise  is 
purely  oscillatory  in  this  situation,  without  attenu- 
ation in  time. 

Equation  (4.15)  or  (4.23)  may  be  integrated  in 
many  other  situations.  A way  of  generating  such 
example  is  discussed  below. 

D.  Properties  of  Vq  (I.Io  I 

A number  of  general  features  exhibited  by 
are  now  observed.  From  (4.15)  we  have 


(4.32) 


tit 


dt 


(4.33) 


which  is  a complex  equation.  The  real  part  of  this 
equation  is  implied  by  (3.2),  but  the  imaginary  part 


I 
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is  an  additional  constraint  on  the  pair  ix(/),  >'(/)  ob- 
tainable from  an  Hq.  Thus  not  every  pair  of  func- 
tions *i(/),  ^'(/)  obeying  (3^)  can  occur  in  a 1/^(1, Iq)- 
Only  those  obeying  (4.33)  are  permitted  in  the  form 
(4.12).  However,  this  does  not  impose  any  con- 
straint on  the  mean- square  fluctuation  behavior  ob- 
tainable in  (3.61),  since  only  |p|,  |v|  are  involved 
here.  Moreover,  (4 .27)  shows  that  any  particular 
pair  of  values  obeying  (3.2)  can  be  obtained 
at  a specific  t with  a proper  choice  of  /„  even  un- 
der (4.3)  with  /,  fixed. 

Any  desirable  pair  of  functions  p(f),  v{t)  obeying 

(4.33)  and  (3.2)  can  be  generated  by  an/,  through 
(4.15).  This  corresponds  to  the  solution  of  the 
Riccati  equation  (4.23)  for  various  /.  For  example, 
we  can  assume  (/,  * w.  fo 

p(f)  cosh(yf),  v(f)  = c***‘’8inh(y/),  (4.34) 

and  determine  the  phase  so  that  (4.33)  is  sat- 
isfied. Equation (4.34)  has  the  interesting  charac- 
teristic that 

<Aa}(f))  = ie^'‘,  (AaJ(0)  = U’'‘,  (4.35) 

without  any  shift  of  phase  to  a'.  Substitution  of 

(4.34)  into  (4.33)  shows  that  0(1)  obeys  the  differen- 
tial equation 

^ 1 

dt  ' co8h’(y<)  + sinh*(y<)  ’ 

so  that 

♦(f)«(«/y)tan*‘[tanh(yO).  (4.36) 

The  corresponding  /,(/)  is  given  by 

/,  = itanh(2yl)- »y/2w.  (4.37) 

Any  arbitrary  state  | ) will  evolve  into 

N,>  = i^o(/.'o)ko>  = / 

(4.38) 

under  Ug.  If  |i/'o>=  l/5o;Po.''o)i  (4.38)  can  be  read- 
ily determined  from  (4.12), 

- k+£(<,<o);»‘(f,/.),i'(f,g)e'»'‘-‘o',  (4.39) 

♦«./,)«  e(i) + 4 • [^*t(i , fo)  1 . (<  •■*0) 


Therefore,  under  an  arbitrary  Ug(f,<o)  a TCS 
will  remain  a TCS  for  all  time,  but  with  p,  p,  and 
v time  varying.  By  solving  p.(f),  i'(/),  and  £(/)  in 
terms  of  Po  ^ •'a  and  p(/,/o).  •'('.O.  and  6(1, la) > 
(4.39)-(4.42)  and  (4.27)  show  explicitly  that  any 
\i3\ii,v)  can  be  obtained  from  any  given 
by  a Ug(t,ta). 

It  also  follows  from  (4.41)  that  if  /,  is  turned  oil 
at  time  /,  then 

p(f)  = p(f, )«'»'*•'/',  v(f)>=v(<,)c'»"-‘/*.  (4.43) 

1 '.us  apart  from  the  optical  oscillation  the  values 
p and  V are  invariant  in  the  absence  of  two- photon 
interaction  in  Hg.  Furthermore,  such  invariance 
is  obtained  when  and  only  when 

/.'O,  (4.44) 

i.e.,  t'(f)°0.  In  general,  p(f)  and  i/(l)  are  deter- 
mined only  by  /,(!)  independently  of  /,(!),  as  seen 
from  (4.15).  This  is  in  accordance  with  the  fact 
that  the  new  quantum  noise  behavior  is  entirely  a 
consequence  of  /,  *0.  Such  correlation  between  /, 
and  n is  also  observed  in  Sec.  IV A. 

On  the  other  hand,  the  added  value  ((1)  depends 
on  both  /,  and  /,  but  is  always  zero  for  /,  ^ 0 , 

i'o(<,<o)|l5)«=|^>,.  when/,  = 0.  (4.45) 

Therefore  all  at  the  f/^’s  are  generated  by  ffg  with 
/,  «0.  The  addition  of  /,  merely  introduces  a con- 
stant shift  to  6,  as  in  (4.5).  The  states  (4.39)  are 
eigenstates  of 

l»(f)»p{f)a(fo)  + i'(f)a»(g,  (4.46) 

where  a(fg)  is  the  Schrodinger-picture  fixed  photon 
annihilation  operator.  We  have 

b(l)  = Ug(t , la)  ada)  U^d , /„)  + £(/)  = aUl . (4 .47) 

Equation  (4.47)  also  shows 
p(l)  aj,t)  + m a^f ) = ado)  + td) 
for  the  Heisenberg  operator 

a,d)»t/J,d,lo)ado)i^od.fo)- 

Equation  (4.48)  implies,  by  solving  for  a,J(t), 
a,d)  = p*d)ado)-  v(f)a»do)-C'd).  (4.49) 

t'd)»vd)t*(/)-P*d)£d).  (4.50) 

Equations  (4.49)  and  (4.39)  give 


where 


v%(t.  /o)  I A,;  Mo.  •'o  > = I A, + £rd); 


A 
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where  n,  v,  £ , and  O are  still  the  solution  to 
(4.1S)-(4.20).  The  state  t/ ^ / , /„) | ;j >,  is  obtained 
from  evolution  of  the  system  baclcwards  in  time 
under  the  same  /{(/)  in  Hg,  and  the  operator  b(i)  is 
similarly  the  ag(/)  moving  baclcwards,  apart  from 
a constant. 

Even  though  no  solution  of  Ug{t,t^)  in  quadrature 
can  be  obtained  in  general,  its  structure  and  be- 
havior have  already  been  fully  determined.  These 
behaviors  are  important  in  the  applications  of 
TCS’s. 

V.  STATE  OF  TWO-PHOTON  EASING  RADIATION 

We  suggest  that  the  output  radiation  of  an  ideal 
monochromatic  two-photon  laser  is  in  a TCS.  In 
contrast,  the  output  of  an  ideal  one-photon  laser  is 
in  an  ordinary  coherent  state.  Thus  TCS’s  are  the 
characteristic  states  of  two-photon  lasing  pro- 
cesses (two  photons  from  the  same  mode). 

A.  Two^hoton  laser 

Two  photons  of  the  same  frequency  w from  the 
same  radiation  mode  can  be  absorbed  in  a single 
atomic  transition  between  two  levels  2 and  1 via 
an  intermediate  state,  i.e.,  a transition  second  or- 
der in  p'  A.  Under  the  usual  dipole  approximation 
the  parities  of  levels  2 and  1 have  to  be  identical. 
First-order  dipole  transitions  between  the  two  lev- 
els are  forbidden  in  this  circumstance.  On  the 
other  hand,  higher-order  multipole  contributions 
to  second-order  p-  A transitions,  which  may  con- 
nect states  of  different  parity,  can  be  significant 
First-order  transitions  between  states  of  opposite 
parity  are  also  possible  via  the  term  in  the 
radiation  interaction  Hamiltonian.^"'^  This  is  pos- 
sible because  in  the  multipole  expansion  of  the 
vector  potential  A(r),  higher-order  multipole  terms 
introduce  atomic  coordinates  into  the  A-  Hamilton- 
ian. Two-photon  absorption  spectroscopy  is  an 
important  experimental  technique  by  now,'’  and 
many  discussions  of  two-photon  absorption  pro- 
cesses can  be  found  in  the  literature. 

The  radiation-matter  (to  be  referred  to  as  atoms) 
interaction  responsible  for  the  two-photon  transi- 
tion can  be  expressed  in  the  form 

H,  = «(pMa^+p*AfV),  (5.1) 

where  p is  a coupling  coefficient  and  M is  the 
atomic  polarization  operator  which  flips  the  state 
of  the  atom  If  the  A^  mechanism  is  under  consid- 
eration, Hi  is  a true  Hamiltonian.'^*"  If  virtual 
transition  through  an  intermediate  state  has  been 
suppressed  through  a nonlinear  susceptibility,  or 
equivalently  through  the  coupling  coefficient  p, 

(5.1)  is  only  an  effective  Hamiltonian.  In  this  case 
it  can  be  derived  by  following  a procedure  similar 


to  the  standard  quantum  theory  of  parametric  in- 
teraction."’''' By  methods  widely  used  in  nonlinear 
optics,'  classical  equations  of  motion  for  the  field 
and  the  active  atomic  variables  can  be  obtained  in 
which  the  virtual  intermediate  states  are  accounted 
for  through  p.  Examples  of  such  treatment  for  two- 
photon  systems  can  be  found  in  many  places.'* 

These  classical  equations  of  motion  with  a quan- 
tum-mechanically  computed  susceptibility  can  be 
directly  quantized.  The  effective  quantum  Hamil- 
tonian follows  from  the  resulting  quantum  equa- 
tions of  motion. 

In  addition  to  absorption,  two-photc’’.  emission 
which  provides  a two-photon  lasing  mechanism  is 
also  described  by  (5.1).  In  considering  the  condi- 
tions for  stimulated  two-photon  emission  in  cither 
the  standing-wave  or  traveling-wave  configuration 
with  loss  included,  it  is  easily  found  that  stimulated 
emission  cannot  be  set  up  from  spontaneous  emis- 
sion alone.  A sufficiently  strong  field  at  the  lasing 
frequency  w has  to  be  present  initially  There- 
fore we  consider  the  two-photon  lasing  configura- 
tions depicted  in  Fig.  2.  It  is  not  the  purpose  of 
this  article  to  present  the  semiclassical  theory  of 
two- photon  laser  operation,  or  to  discuss  the  op- 
erating consideration  in  a realistic  material  sys- 
tem. But  the  following  facts  may  be  observed: 

In  the  traveling- wave  configuration.  Fig.  2(a), 
the  gain  is  obtained  from  a large  external  pump 
field  which  may  have  many  frequencies  but  which 
is  still  in  a single  temporal  mode.  Neglecting  var- 
ious loss  mechanisms,  this  lasing  process  can  be 
described  by  the  interaction  Hamiltonian  (5.1)  with 
an  initial  radiation  state  having  a nonzero 
amount  of  power,  ttp^a^a  eO.  If  a laser  is  used  for 
the  external  pump,  = | a ) (a  [ , but  intense  ther- 
mal light  can  also  be  used.  The  radiation  state  p, 
as  well  as  other  dynamical  variables , are  spatial- 
ly dependent  in  this  case.  In  the  standing-wave 
case.  Fig.  2(b),  only  a single  standing-wave  mode 
is  under  consideration.  The  gain  is  initially  ob- 
tained from  the  usual  one-photon  lasing  mechanism, 
the  host  atoms  for  both  lasing  transitions  being 
contained  in  the  same  cavity.**  The  Hamiltonian  is 
now 


(0)  Ej-E|  • 2li»  (61  Ej-E|  ; Jti. 

Ei-E,'  ■ ti. 


FIG.  2.  Two  possible  configurations  for  stimulated 
two-photon  emission. 
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H = K(u>a^a + +/>,  A/  Ja  ’ + -^p^M  [o) , 

(5^) 


through  the  Fokker-Planck  equations  for  p‘"*(a*,  a). 
In  parttcular,  all  of  the  characteristic  functions 
are  still  in  Gaussian  form,*"  and  for  y*4r, 


with  a vacuum  initial  state;  the  subscripts  2 and  1 
for  p and  M refer  to  the  two-photon  and  one-photon 
lasing  material.  Population  inversion  is  required 
in  both  configurations.  Other  possible  configura- 
tions for  two-photon  lasers  can  be  developed,  but 
these  two  seem  to  be  the  most  important . 

Spontaneous  two-photon  emission  has  been  ex- 
perimentally observed but  not  stimulated  two- 
photon  emission.  Two-photon  stimulated  emission 
was  first  suggested  for  giant  pulse  generation.'*’ 

When  such  a two-photon  laser  operates  far  above 
threshold  with  amplitude  and  phase  stability,  the 
fluctuation  in  M should  be  relatively  small,  so  that 
M is  approximately  a c number.  By  ignoring  the 
quantum  nature  of  Af,  (5.2)  or  (5.3)  becomes  an  Hq 
of  the  form  (4.2).  If  the  dissipative  mechanisms 
are  also  neglected,  (4.12)  shows  that  a TSC  |/3), 
will  be  obtained.  A similar  argument  was  used  by 
Glauber**  to  suggest  that  a coherent  state  describes 
one- photon  laser  radiation  far  above  threshold.  As 
we  show  in  the  Appendix,  this  can  be  formalized  in 
the  so-called  self-consistent-field  approximation 
(SCFA),  which  has  been  used  in  one- photon  laser 
theory .”•**  Therefore  a two-photon  device  of  this 
type  operating  well  above  threshold  can  be  expected 
to  produce  a state  with  further  small  classi- 
cally random  fluctuations  in  the  varL''.''!c  P,  Ln  a 
way  exactly  analogous  to  one-photon  lasers  where 
I a ) is  produced  with  small  classical  fluctuations 
in  a.  This  argument  can  be  made  for  a standing- 
wave,  a traveling- wave,  or  an  arbitrary  spatial- 
temporal  mode. 


B.  Effect  of  perturbations  in  the  two-photon  system 

An  additive  Gaussian  noise  and  a damping  me- 
chanism can  be  introduced  to  account  for  some  un- 
avoidable sources  of  noise  similar  to  the  one-pho- 
ton case,  while  still  preserving  that  distinguishing 
feature  of  TCS’s  with  and<Aa'*)  smaller  than 

This  can  be  illustrated  as  follows:  The  quan- 
tum effect  of  loss  is  incorporated  through  the  den- 
sity operator  relaxation  equation 

I?  *-T[^fo,p]+aY(2apa^-a’ap-paM 
9t  h 

+ Yfi{a*pa  + apd*  - a^ap- paa^),  (5.3) 

with  Hq  given  by  (4,2).  The  damping  rate  is  y and 
Jf  represents  the  thermal  noise  contribution  from 
the  reservoir  that  couples  to  the  radiation  field. 
This  quantum  description  of  loss  has  been  exten- 
sively discussed  in  the  literature. For  an  ini- 
tial coherent  .state,  (5.3)  can  be  solved  exactly 


(Aa«,(f)>=l 


Y 

yi4r 


r 

^y±4r 


yff  It  --lrnr)t\ 

+ ~ ® It 

2(y±Ar) 


(5,4) 


with  a'  = ae'“*.  It  can  be  readily  shown  that  <Aa5*(f)) 
>0  and  that  the  uncertainty  principle  is  obeyed  for 
all  1.  Note  that  even  when  R = 0,  the  state  is  no 
longer  pure  in  the  presence  of  y.  Asymptotically 
we  have 


<Ao,'»(/))-r>f/8*'.  <Aa,'*(f))-«';  t—o,  2r»y. 

(5.5) 


The  condition  2r»r  can  be  satisfied  together  with 
the  stability  condition  (4.3)  for  the  usual  values  of 
y. 

u3>2r»y.  (5.6) 


Thus  for  If  ^ 1 the  corresponding  state  approaches 
an  eigenstate  of  a, . In  any  case  the  quantum  fluctu- 
ation in  a^  vanishes  as  < — «» . 

From  the  results  of  Sec.  IV  B,  it  can  be  seen  that 
two-photon  laser  is  a phase- sensitive  device.  If  a 
uniform  random  phase  is  assumed,  i.e.,  if  a'  = ae** 
for  <p  uniformly  distributed  in  [-*,  vj  is  the  dy- 
namical observable  one  has  to  deal  with,  then 

<Aa«>  = <Aa«>  = i«AaJ>  + <Aa*».  (5.7) 

Equation  (5.7)  contributes  a noise  greater  than  J 
lor  both  quadrature  components  and  one  distin- 
guishing characteristic  of  a TCS  is  lost.  The  de- 
tailed theories  of  one-photon  lasers  show  that  am- 
plitude fluctuation  is  suppressed  in  a laser  operat- 
ing far  above  threshold , and  the  phase  fluctuation 
is  also  inversely  proportional  to  the  average  pho- 
ton numbers.  Moreover,  these  amplitude  stabili- 
zation and  phase  stabilization  properties  are  gen- 
eral features  of  a wide  class  of  self- sustained  os- 
cillators, quantum  or  classical.”*’®  If  the  two- 
photon  laser  belongs  to  this  class,  the  phase  would 
diffuse  away  from  some  initially  measured  value 
only  slowly  for  high-power  opeiation  similar  to  the 
one- photon  case.  The  phase  can  then  be  tracked 
and  the  relative  phase  determined  accurately.  As 
a result  the  field  may  be  regarded  to  be  in  a true 
TCS.  Whether  the  perturbation  introduced  by  the 
atomic  fluctuations  may  indeed  be  sufficiently 
small  so  that  | /3  >,  does  not  become  | p ) can  be  de- 
termined only  by  explicitly  working  out  the  detailed 
quantum  theory  of  two-photon  lasers.  The  quantum 
theory  of  one-photon  lasers  has  been  extensively 
developed  but  the  two- photon  laser  requires 

more  careful  approximations  because  of  the  small 
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noise  in  one  quadrature  component  and  the  more  itial  states  of  the  atoms  are  identical,  but  they  are 

delicate  equations  of  motion  involved.  Hopefully,  different  in  the  two-photqn  case.  The  question  of 

the  TCS  picture  can  be  checked  experimentally,  as  quantum  statistics  is  similar  in  both  cases.  The 
discussed  in  Sec.  VI,  by  practically  realizing  a pump  field  of  Fig.  3 should  be  treated  as  a quan- 

two- photon  laser.  In  any  case,  | P ),  serves  as  the  turn  degree  of  freedom.  A fully  quantum  treatment 

basic  description  to  be  averaged  further  by  a of  degenerate  parametric  oscillators  can  be  equal- 

classical  density  function  in  the  presence  of  other  ly  or  more  complex  than  a fully  quantum  treatment 

random  interference.  of  two-photon  lasers,  depending  on  whether  the 

Recently,  the  quantum  theory  of  multiphoton  atomic  variables  are  included  in  the  parametric 

lasers,  including  a quantum  treatment  of  the  atoms  case  and  whether  the  important  intermediate  state 

and  the  loss  n.echanism,  was  discussed  by  McNeil  is  included  in  the  two-photon  case  when  present, 

and  Walls.”"®*  However,  the  nature  of  their  ap-  No  single-mode  degenerate  parametric  oscilla- 

proximations  was  not  clarified  and  they  produced  tion  above  threshold  has  been  observed.  Careful 

lasing  results  from  models”’®  which  cannot  lase  consideration  has  to  be  given  to  any  possible  com- 

because  of  the  absence  of  an  input  field  to  provide  peting  processes  in  the  actual  development  of  a 

a nonzero  gain.  While  some  qualitative  features  two-photon  laser  or  degenerate  parametric  ampll- 

of  their  two-photon  results  are  in  agreement  with  fier,  such  as  third- harmonic  generation, 
the  TCS  picture,  such  as  large  total  fluctuation 

from  (3.35),  they  suggest  that  the  photon  statistics  yi  COHERENCE  AND  PHOTON  COUNTING 

are  similar  to  the  one-photon  case.  The  Poisson- 

like  distribution  they  find  in  Ref.  59  is  obtained  for  Two-photon  coherent  states  are  not  coherent  ac- 

an  unrealistic  loss  model,  in  addition  to  the  above  cording  to  Glauber’s  original  definition.”  That  is, 

error.  The  lossless  photon  distribution  they  ob-  they  do  not  factorize  the  correlation  functions 

tain  in  Ref.  61  is  neither  Poisson  nor  of  the  form  which  describe  photon  counting  experiments  em- 

(3.23).  Since  the  atomic  populations  rather  than  ploying  ordinary  photodetectors.  However,  they 

the  atomic  polarizations  were  treated  as  a con-  factorize  a different  set  of  correlation  functions 

stant  in  the  derivation  of  that  distribution,  it  ap-  obtained  by  transforming  a to  6 in  the  field  expan- 
pears that  it  is  valid  near,  instead  of  far  above,  sion.  This  kind  of  transformed  correlation  func- 

threshold.  A complete  and  careful  calculation  of  tion  gives  the  joint  counting  rates  for  detectors 

the  two- photon  laser  state  remains  to  be  carried  operating  with  an  interaction  mechanism  different 

from  the  usual  dipole  absorption,*"  in  analogy  with 
detectors  that  function  by  stimulated  emission 

discussed  by  Mandel.®"  Thus  a TCS  is  fully  co- 

C.  Degenente  primetric  .mplifier  ^ 

A degenerate  parametric  amplifier  interpretation  finition  of  optical  coherence.  Further  detailed  dis- 
can be  given  to  (4.2)  similar  to  the  above  two-pho-  cussion  can  be  given  only  in  the  multimode  t-eat- 

ton  laser  interpretation.  When  the  material  sys-  ment  of  TCS. 

tern  is  suppressed  through  a susceptibility  and  Here  let  us  consider  photon  counting  experiments 

when  the  pump  field  is  treated  classically,  a two-  associated  with  a single  radiation  mode  whose 

photon  coherent  state  would  be  obtained  from  spatial  dependence  has  been  integrated  through  the 

a degenerate  parametric  amplifier  such  as  the  one  spatial  response  of  a photoabsorption  detector . 

depicted  in  Fig.  3.  An  intense  initial  field  at  Ire-  For  any  state  p the  distribution  <n|p|n>  can  be 

quency  w has  to  be  present  in  addition  to  the  pump  generated  by  an  ideal  detector  if  the  counting  time 

field  for  degenerate  parametric  amplification.  A F is  large  enough.®’  Therefore  we  consider  the 

TCS  is  already  present  in  the  treatment  of  nonde-  generating  function 

generate  parametric  amplifiers  by  Mollow  and 

Glauber"’  in  the  joint  Wigner  distribution  for  the  Q(x)=^(l  - X)''<n|p|«>,  (6.1) 

signal  and  idler  modes,  even  though  a coherent  ^ ’ 

state  plus  noise  is  obtained  for  each  of  these 

modes.”  That  a TCS  can  be  generated  in  a single-  = (N- • • (N- !)•  • • (JV - n + 1)>  (6.2) 

mode  degenerate  parametric  amplifier  was  also  i 

observed  by  Stoler”  and  by  Lu.*®  8X"  j ' 

The  phase- sensitive  behavior  of  a two-photon 

laser  is  reminiscent  of  parametric  processes,  and  With  the  Glauber  model'®  of  an  n„-atom  photodetec- 

a degenerate  parametric  amplifier  is  a kind  of  two-  tor  with  quantum  efficiency  s and  a sufficiently 

photon  laser.  A basic  distinction,  however,  can  small  counting  time  T , the  generating  function  for 

be  made.  In  the  parametric  case  the  final  and  in-  the  counting  distribution  l>{n,T)  of  registering  n 
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counts  In  time  T is  merely  a scaling  of  (6.1), 


(-1)"  9"Q(sw,.Tx) 
n!  ax" 


(6.4) 


From  (3J23)  the  function  Q(x)  for  a TCS  can  be 
evaluated  in  closed  form  using  the  summation  for- 
mula for  Hermite  polynomials,’’ 


Q(X)»T*''‘exp^|.(l-  X)t-  1]  1p1’ 


(3 


(6.5a) 

t.[|m|‘-(1-x)*|i'1>1-‘.  (6.5b) 

From  (6.5)  the  coefficients  = (a^"a">/<oM"  a 
first-order  coherent  field’*  can  be  computed  by 
(6.3).  In  particular,  <a’a>=  |«'|*+ 14^1’,  as  given  by 
(3.29),  and 

(a«a*>»<a»a>’  + 2i^l’(4|v|»  + 3)|i/|’ 

- (p*i'*^  + pv6*’)(l+4lv|’)+  |v|’  + 2lvl". 

(6.6) 

Any  value  >0  can  be  generated  from  (6.6)  by 
proper  choice  of  (p,t',^).  In  particular,  <1  for 
a TCS  if  and  only  if 

2|i/|«+|i.|*+2(4|v|‘  + 3)|i/|»|/3|’ 

<(p*v*0*+pi»0*»)(l  + 4|v|’),  (6.7) 


and  1 if  the  sign  is  reversed.  For  0 = 0,  we 
always  have^j>l  from  the  noise  ji'l’,  but^j<l  is 
clearly  possible  for  small  v and  large  p with  non- 
vanishing 0.  Recall  that  ^2  < 1 is  a characteristic 
quantum  effect  without  classical  analog,  and  that 
currently  available  sources  all  generate g^^l. 
Stoler”  and  Lu‘®  have  previously  observed  that  gj 
<1  may  be  obtained  from  a state  of  this  kind. 

For  most  values  of  (PiV),  g^  can  be  tuned  by 
varying  0 in  a TCS  from  (6.7),  exhibiting  both  anti- 
correlation (^2'^1)  enhanced  correlation  (gj 
>2)  effects.  This  fact  permits  one  to  check  whether 
a certain  particular  TCS  is  obtained  by  a Hanbury- 
Brown-Twiss  experiment.  The  photon  counting 
distribution  (3.23)  or  (6.4)  is  also  sufficiently  dif- 
ferent from  a Poisson  distribution  that  its  measure- 
ment will  reveal  the  presence  of  a TCS.  The  in- 


u 


FIG.  3.  Phase-matched  degenerate  parametric  amp- 
lification with  pump  field  at  frequency  2u. 


fluence  of  nonideal  detectors  and  noise  on  these 
photocount  statistics  should  be  worked  out  in  de- 
tail for  actual  comparison  with  experiments. 

VII.  APPLICATIONS  OF  TWO-PHOTON 
COHERENT  STATES 

We  will  briefly  describe  some  novel  applications 
of  two-photon  lasing  processes  which  result  from 
the  quantum  properties  of  TCS.  Quite  apart  from 
rheir  quantum  characteristics,  two-photon  lasers 
also  exhibit  certain  pulse- shortening  behavior  ab- 
sent in  the  one-photon  case.”  Their  particular 
nonlinearity  may  lead  to  many  interesting  phenom- 
ena and'applications,  including  the  realization  of 
very-high-power  lasers.  However,  only  certain 
applications  of  TCS’s  will  be  treated  here.  These 
applications  also  illustrate  the  essential  differences 
between  TCS’s  and  ordinary  coherent  states. 

A.  Ideal  linear  amplifier 

Under  the  action  of  a two-photon  amplifier,  an 
input  coherent  state  |0)  will  be  changed  to  |0)2, 
from  (4.45).  For  simplicity  we  assume  p and  v to 
be  real.  From  (3.28), 

(ai>r'0l'={M- «')0i,  (7.1a) 

We  ”4 '<!’  + •')  0,.  (7.1b) 

Thus  one  of  the  original  quadrature  components  is 
linearly  amplified  while  the  other  is  linearly  at- 
tenuated. The  nonlinearity  of  a two- photon  ampli- 
fier is  therefore  of  the  simplest  kind;  it  is  piece- 
wise  linear. 

To  be  definite  let  v>Q.  The  quantum  noise  (A0J) 
decreases  from  (A0J)  = j to  (A0j)  = j(p-  ix)’,  while 
that  of  02  increases  to  ?(m  + i')“.  from  (3.31).  If  we 
define  the  signal-to-noise  ratio 

(S/AV  = (a//(Aa5>,  ; = 1,2,  (7.2) 

where  the  average  is  taken  in  a state  p,  we  find 
that  (S/N)j  is  invariant  in  the  amplification  and  at- 
tenuation, i.e., 

{S/JV)p„  = (S/N)„„^.  (7.3) 

The  above  definition  of  signal-to-noise  ratio  is 
meaningful.  If  we  try  to  determine  the  value  of 
by  measuring  02  in  s’  coherent  state  or  TCS,  (7.2) 
provides  a standard  or  measure  on  how  well  one 
may  estimate  02,”’*  or  equivalently  how  "noisy”  0, 
is  in  the  "best”  measurement. 

We  cw  interpret  |0)  as  the  state  of  an  optical 
signal  presented  to  a receiver  which  generates 
|0>,  by  a two-photon  amplifier  and  then 

measures  a,  o"  |(3)f  H is  obvious  physically  that 
the  signal-to-noisc  ratio  cannot  be  Improved  by 
amplification.  The  introduction  of  into  t'g  does 
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not  increase  the  signal  energy,  i.e.,  the  part  of 
which  is  proportional  to  The  advantage  of  am- 
plification is  practical,  similar  to  the  classical 
case.  Raising  the  signal  level  makes  possible 
easier  measurement  of  (Jj.  Since  the  two-photon 
amplification  process  (4.45)  operates  above  thres- 
hold, it  can  be  superior  to  one-photon  amplifiers 
operating  below  threshold  in  its  possibilities  of 
yielding  higher  gain  and  lower  noise.  If  amplifica- 
tion for  both  quadrature  components  is  desired, 
the  input  beam  can  be  first  split  and  two  different 
VgUtO  applied  to  di  and  *1,  separately.  Ideally, 
no  signal-to- noise  degradation  will  occur  from 
(7.3). 

The  above  amplifier  performance  is  not  sensitive 
to  the  addition  of  a Gaussian  noise.  Since  we  are 
now  utilizing  the  bigger  noise  component,  there  is 
no  need  to  maintain  <AaJ)<{  and  a random  phase 
would  not  seriously  degrade  the  system  perfor- 
mance. Therefore  the  process  described  by  f/o(/,/o) 
provides  a possible  mechanism  in  the  realization 
of  an  ideal  linear  amplifier. 

B.  Local  oscillator 

One  may  take  advantage  of  the  absolutely  small 
<AaJ)«i  in  \ fi),  in  the  reception  of  a coherent- 
state  signal  with  a size-limited  detector.  By  em- 
ploying a local  oscillator  which  generates  radiation 
In  state  it  is  possible  to  attenuate  the  quan- 
tum noise  by  an  amount  which  compensates  the 
Urge  diffraction  loss  incurred  in  free  space  or  un- 
guided  propagation.  This  is  achieved  by  coupling  a 
mode  in  to  the  incoming  mode  in  |a)  in  the 
form  of  a new  mode 

a,  = €as  + (l-<»)‘'X.  (''•4) 

where  c* represents  the  fraction  of  energy  received 
and  ttg,  and  og  are  the  photon  annihiUtion  op- 
erators for  the  signal,  the  local  oscillator,  and 
the  receiver  mode.  The  quantum  noise  in  Og  is 
clearly  dominated  by  that  of  for  small  € The 
signal-to-noise’  ratio  in  og,  becomes  ~i  a J for 
(Aa^j>=«0,  whereas  it  is  for  a coherent- 

state  local  mode  <Anii)  = h jTjtfttBr  discussion 
can  be  found  in  Ref.  1.  Na^MPt^KIS  f^yantage  of 
i cannot  be  obtained^lf  a one-photon- laser 
local  oscilUtor  is  used  to  produce  a cohwent- 
state  receiver  mode  Og  and  a f/g(/,/o)  subsequently 
applied,  because  of  the  corresponding  signal  at- 
tenuation. An  absolutely  smaU  is  needed  in 

this  scheme. 

If  an  arbitrary  amount  of  power  is  available,  it 
is  possible  to  have  a state  with  (Aaj)»0. 
ever,  for  a fixed  total  radiation  energy 

trprt^a^S,  (7.5) 


a decrease  of  (Aaj)  in  Id),  can  be  obtained  only  at 
the  expense  of  spending  a portion  of  available  en- 
ergy S in  the  form  of  added  quantum  noise  energy 
lv|*,  from  (3.35).  Thus  (Aaj)-O  requires  S-”. 
This  consideration  is  not  important  in  the  context 
of  local  oscillators  at  a receiver  where  large 
enough  amount  of  power  is  usually  available  to 
make  (Aa’)  sufficiently  small.  A four-orders-of- 
magnitude  reduction  from  the  quantum  noise  i of 
a,  (Ivies' 10^)  for  a GHz  bandwidth  optical  signal 
at  u/2)t=^10‘’  Hz  requires  only  a Ivj-  correspond- 
ing to  ~1  pW.  This  will  usually  bring  it  down  to  the 
level  of  other  extraneous  noises.  Furthermore, 
even  a ten-orders-of-magnitude  reduction  of  the 
qpiantum  noise  requires  -1 W under  the  same  con- 
dition. 

Other  practical  considerations  also  should  not 
invalidate  the  utilization  of  this  small-noise  prop- 
erty. With  proper  signaling  and  reception  schemes 
one  can  always  observe  the  signal  quadrature  com- 
ponent without  significant  degradation  from  ran- 
dom-phase modulation.  From  (4.43)  it  is  assured 
that  free-space  propagation  of  1*3)^  will  not  disturb 
the  quantum  noise  characteristics.  The  rotation 
(3 .59)  can  be  readily  compensated  in  an  optical 
receiver  without  affecting  the  signal  level,  so  that 
TCS’s  are  as  good  as  minimum-uncertainty  states 
in  this  application.  These  and  many  other  problems 
involved  in  a realistic  communication  system  will 
be  treated  elsewhere. 

C.  Transmission  source 

Further  performance  improvement  can  be  ob- 
tained if  the  transmitter  generates  information- 
carrying radiation  in  a TCS.  In  this  case  the  radi- 
ation power  constraint  (7.5)  poses  a major  limit- 
ation on  the  achievable  signal-to-noise  ratio.  It 
has  been  shown  that®’  TCS's  provide  the  maximum 
{S/N)i  obtainable  for  one  of  the  quadrature  compo- 
nents among  all  possible  states  satisfying  (7.5). 

By  spending  a traction  S/(2S  + 1)  of  S as  quantum 
noise  energy  with,  e.g., 

p=(S  + l)/(2S  + l)-•'^  »'-S/(2S  + l)-‘/»,  (7.6) 

Ps  = 4s  = 0,  (7.7) 

we  have 

(S/^)ii.>,  = 4(S*  + S).  ‘ (7.8) 

Compared  to  the  signal-to-noise  ratio 

(SAV),,g,  = 4S  (7.9) 

in  a coherent  state,  we  see  that  the  available  "sig- 
nal energy”  is  effectively  increased  from  S to  S- 
+ S,  This  is  actually  a reduction  of  quantum  noise.®* 
Among  other  things,  this  increase  in  (S/.V),  leads 
to  a higher  information  capacity  even  when  the 
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other  quadrature  component  is  not  employed  for 
information  transmission.*  In  the  presence  of 
other  extraneous  noise,  quantum  noise  reduction 
in  a transmitter  state  |(i),  can  also  be  obtained 
with  only  a relatively  small  Iv]’,  as  illustrated 
above. 

Finally,  we  observe  that  the  new  photocount  sta- 
tistics may  itself  lead  to  novel  applications  in 
communications. 

VIII.  CONCLUDING  REMARKS 

We  have  discussed  in  some  detail  the  properties 
of  a class  of  radiation  states  which  are  closely 
related  to  the  well-known  coherent  states.  These 
two-photon  coherent  states  possess  a number  of 
distinguishing  characteristics  which  lead  to  poten- 
tial device  applications  for  quantum  noise  reduc- 
tion. They  may  also  find  application  in  experi- 
mental situations  that  require  low  noise  sensitivity. 

It  has  also  been  suggested  that  these  states  can 
be  obtained  from  two-photon  stimulated  processes, 
or  equivalently  that  they  provide  a good  description 
of  two-photon  laser  states.  The  ultimate  validity 
of  this  proposition  must  be  tested  by  an  elaborated 
quantum- statistical  theory  of  two-photon  lasing 
processes  and  above  all,  hopefully,  by  experiment. 

Our  present  work  was  motivated  by  the  desire  to 
realize  physically  certain  specific  quantum  states 
and  quantum  measurements  for  practical  applica- 
tions. It  appears  that  the  general  problems  of 
state  and  measurement  synthesis  are  interesting 
and  useful  areas  in  quantum  optics  and  electronics, 
as  well  as  in  the  foundations  of  quantum  mechanics. 
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APPENDIX:  TWO-PHOTON  LASING  STATE 
IN  THE  SCFA 

Here  it  is  shown  that  the  radiation  state  of  a two- 
photon  laser  is  a TCS  in  the  self-consistent-field 
approximation  (SCFA),”’**  if  loss  is  also  neglected. 
Let  H be  the  total  Hamiltonian  of  an  atom-field 
system; 

(Al) 

with  given  by  (5.2).  In  the  SCFA,  the 

total  density  operator  is  assumed  to  factorize, 

P = Pa»®P/d-  (A2) 

■ The  equation  of  motion  for  the  radiation  density 
operator  can  be  obtained  from  that  of  p by 
tracing  over  the  atomic  variables  with  the  result 

^=-^[^0.  P/dl,  (A3) 

where  /fg  is  given  by  (4.2),  with 

/i  = ‘^.  /2=Patr(A/,p,,),  /j=p,tr(,Vf,p„).  (A4) 

Thus  a TCS  will  be  developed  from  an  initial  co- 
herent state  from  (4.12),  and  (A4)  gives  the  effec- 
tive field  Hamiltonian  under  the  SCFA . Note  that 
a resonant  factor  is  automatically  present  in 

/j,  from  the  atomic  resonance  in  M. 

Even  though  the  above  derivation  is  given  for  a 
single  frequency,  it  applies  equally  to  any  spatial- 
temporal  mode  if  is  in  the  form  (5.1).  The 
same  result  is  obtained  for  a parametric  oscillator 
by  interpreting  as  the  pump  photon  destruction 
operator  and  p^  as  the  susceptibility,  with  Af,  = 0. 
However,  it  seems  less  justified  to  neglect  the 
quantum  fluctuation  of  the  pump  photon,  which  is 
also  described  by  a boson  operator  as  the 
signal  mode. 
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Abstract 

Recent  theoretical  work  has  shown  that  novel  quantum  states, 
called  two-photon  coherent  states  (TCS),  have  significant  poten- 
tial for  improving  free-space  optical  communications.  In  this 
three-part  study,  the  communication  theory  of  TCS  radiation  is 
explored.  The  present  work,  part  I,  addresses  the  Issues  of  quantum- 
field  propagation  and  optimum  quantum  state  generation.  In  particular, 
the  quantum  analog  of  the  classical  Huygens-Fresnel-Klrchhoff 
diffraction  theory  for  quaslmonochromatlc  scalar  waves  is  developed. 
This  result,  which  describes  the  propagation  of  arbitrary  quantum 
states  as  a boundary-value  problem  suitable  for  communication-system 
analysis,  is  used  to  treat  a number  of  quantum  transmitter  optimiza- 
tion problems.  It  is  shown  that  under  near-field  propagation  con- 
ditions a TCS  transmitter  maximizes  field-measurement  slgnal-to- 
noise  ratio  among  all  transmitter  quantum  states;  the  performance 


I 

of  the  TCS  system  exceeds  that  for  a conventional  (coherent 
state)  transmitter  by  a factor  of  N + 1,  where  N Is  the  average 
number  of  signal  photons  (transmitter  energy  constraint).  Under 
far-fleld  propagation  conditions.  It  Is  shown  that  use  of  a TCS 
local-oscillator  In  the  receiver  can.  In  principle,  attenuate 
field-measurement  quantum  noise  by  a factor  equal  to  th<-  diffrac- 
tion loss  of  the  channel.  If  appropriate  spatial  mode-mixing  can 
be  achieved.  These  communication  results  are  derived  assuming  that 
field-quadrature  quantum  measurement  Is  performed.  In  part  II 
of  this  study,  photoemlsslve  reception  of  TCS  radiation  will  be 
considered;  It  will  be  shown  therein  that  homodyne  detection  of 
TCS  fields  can  realize  the  field-quadrature  slgnal-to-nolse  ratio 
performance  of  part  I.  In  part  III,  the  relationships  between 
photoemlsslve  detection  and  general  quantum  measurement  are  estab- 
lished. In  particular,  a synthesis  procedure  will  be  obtained  for 
realizing  all  the  measurements  described  by  arbitrary  TCS. 

I 
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I Introduction 

* ••The  quant urn- me chan leal  theory  of  optical  communications  • 
seeks  to  delineate  the  ultimate  physical  constraints  Imposed 
on  the  transmission  of  Information.  The  work  to  date  In  this 
field  has  focused  on  determining  optimum  quantum  measurements  for 
various  detection  and  estimation  problems  [l]-[4].  However, 
to  determine  the  ultimate  performance  limitations  Imposed  by 
quantum  effects.  It  Is  also  essential  to  consider  optimum  quantum- 
state  generation,  because  post-measurement  statistics  (hence 
detection  and  estimation  performance)  depend  on  both  the 
quantum  state  and  the  quantum  measurement.  Indeed,  Yuen  has 
shown  that  certain  ’’generalized”  coherent  states  of  the  radia- 
tion field  possess  novel  quantum  noise  characteristics  which 
offer  the  potential  for  greatly  Improved  optical  communications 
C5]-[7].  These  new  states  are  related  to  the  ordinary  co- 
herent states  of  Glauber,  but  have  not  been  experimentally 
observed.  They  have  been  called  two-photon  coherent  states, 
because  they  can  be  generated.  In  principle,  by  stimulated 
two-photon  processes  [?]•  In  this  paper  and  the  sequel  [8],  [9], 
we  shall  consider  the  use  of  two-photon  coherent  state  (TCS) 
radiation  In  free-space  optical  communications. 

In  the  usual  formulation  of  quantum  communication  theory, 
the  electromagnetic  field  Is  quantized  by  representing  the 
- positive-frequency  field  within  a receiver  region  as  a sum  of 
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orthogonal  space-time  modes,  and  ImposlnR  the  photon  annihila- 
tion operator  commutation  rule  on  the  coefficients  In  this 
expansion  [M,Ch.5],  [10],  This  Initial-value  problem  approach  Is 
not  well-suited  to  the  study  of  optimum  quantum  transmitter 
design.  For  such  a study.  It  Is  necessary  to  have  a quantum- 
mechanical  diffraction  theory,  analogous  to  the  classical  Huygens- 
Presnel-Klrchhoff  theory,  in  which  quantum  field  propagation  In 
vacuum  can  be  rigorously  treated  as  a boundary-value  problem. 

Because  some  of  the  Issues  of  a quantum  diffraction  theory  are  not 
directly  related  to  communication, we  will  not  present  a full 
treatment  of  such  a theory  In  this  paper.  On  the  other  hand,  we 
shall  develop.  In  Section  II,  those  features  of  quantum  state 
propagation  which  provide  the  basis  for  representing  any  optical 
channel  In  a quantum-mechanical  form  convenient  for  communication 
analysis. 

The  development  of  quantum  diffraction  theory  will  be  for 
arbitrary  quantum  states;  in  Sections  III-V  we  shall  apply  this 
theory  to  TC6  radiation.  Under  near-fleld  propagation  conditions, 
with  a fixed  field-quadrature  measurement.  It  will  be  shown  that  a 
TC6  transmitter  maximizes  post-measurement  slgnal-to-nolse  ratio 
over  all  transmitter  quantum  states  satisfying  an  average  energy 
constraint.  The  TCS  system  substantially  outperforms  a conventional 
transmitter  of  the  same  energy.  Under  far-fleld  propagation 
conditions,  '■rlth  a fixed  field-quadrature  measurement.  It  will  be 
shown  that  use  of  a TCS  local-oscillator  In  the  receiver  can.  In 
principle,  attenuate  post-measurement  quantum  noise  by  a factor 
equal  to  the  diffraction  loss  of  the  channel  * If  appropriate  spatial 
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mode-mixing  can  be  achieved. 

In  deriving  the  foregoing  communication  results,  we  shall 
assume  that  the  field-quadrature  quantum  measurement,  i.e.,  that  of 
a^  the  cosine  quadrature  of  the  TCS  radiation  mode,  is  performed. 
When  using  conventional  radiation,  the  a^  measurement  can  be 
accomplished  by  homodyning  [11].  However,  because  the  statistics 
of  photoemissive  detection  depend,  non-trlvlally , on  the  quantum 
state  of  the  radiation  field,  the  preceding  statement  regarding 
homodyning  must  be  re-examined  for  TCS  sources.  In  [8]  we  shall 
present  a detailed  analysis  of  photoemissive  TCS  detection; 
homodyning  will  be  shown  therein  to  realize  the  same  signal-to- 
nolse  ratio  as  the  a^  measurement  for  an  arbitrary  quantum  state. 
The  full  equivalence  of  homodyning  and  the  a^^  measurement  will  be 
established  in  [9]. 
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II  Quantum  Field  Propagation 

Open  almost  any  graduate  text  on  communication  theory  Cl2]“Cl^], 
and  there.  In  the  Introductory  chapter,  you  will  find  a generic 
block  diagram  for  a single  communication  link,  much  like  that 
shown  In  Fig.  1.  The  ubiquity  of  Fig.  1 stems  from  Its  Inclusion 
of  virtually  all  single-link  communication  problems,  e.g.,  source 
coding,  signal-set  selection,  channel  coding,  and  optimum  receiver 
principles.  Its  appearance  In  the  1970  review  article  on  quantum 
communication  theory  by  Helstrom  et  al . [1],  belles  the  fact 
that  to  that  time  attention  had  been  focused  almost  exclusively 
on  receiver  optimization.  Moreover,  Helstrom's  recent  monograph  [4] 
attests  to  the  continued  absence  of  a quantum-transmitter  op- 
timization theory. 

A primary  cause  for  the  present  dearth  of  quantum-transmitter 
analysis  Is,  we  believe,  the  cumbersome  quantization  procedure 
employed  In  prior  studies.  The  radiation  field  Is  usually  quantized 
at  some  Initial  time  t^  within  a spatial  volume  V that  Is  large  enough 
to  contain  the  field  over  the  space-time  region  of  Interest  [10 ]. 

In  order  to  use  the  volume  quantization  procedure  to  develop 
a quantum  transmitter  theory,  we  must  take  V to  be  a very  large 
spatial  volume  containing  the  transmitter,  the  propagation  path, 
and  the  receiver.  Once  the  field  has  been  quantized.^  at  an  appro- 
priate Initial  time,  within  V the  quantum  characteristics  of  the 
receiver  entrance— aperture  field  over  the  detection  Interval  can 
In  principle  be  determined. 


^ . 1 

The  foregoing  volume-quantization  procedure  for  quantum 
field  propagation  [15]  Is  seldom  employed.  When  all  the  radiation 
sources  produce  fields  which  are  In  coherent  states  or  their 
classlcally-random  superpositions,  the  radiation  field  In  the 
I neighborhood  of  the  receiver  can  be  determined  by  classical 

! diffraction  theory  techniques  (c.f.  Section  11,2),  On  the  other 

hand,  for  TCS  radiation  ( and  other  quantum  states  which  do  not 
possess  well-behaved  P-representatlons) , classical  techniques  do 
[ not  suffice.  Moreover,  to  use  the  volume-quantization  procedure 

' In  this  case  we  must  explicitly  describe  the  radiation  source 

In  quantum-mechanical  terms.  The  need  to  Introduce  the  source 
dynamical  variables  Into  the  propagation  problem  unduly  complicates 
the  task  of  explicitly  determining  the  quantum  characteristics 
of  the  received  field.  Thus,  our  first  concern  herein  will  be  to 
formulate  a quantum  diffraction  theory,  similar  In  spirit  to  Its 
classical  counterpart.  In  order  that  the  quantum  field  generation 
and  field  propagation  problems  can  be  separated.  Quantum  channels 
can  then  be  conveniently  analyzed  within  the  framework  of  Pig,  1, 

II, 1 Specification  of  Quantized  Aperture  Fields 

Classical  diffraction  theory  treats  electromagnetic  wave 
propagation  as-  a boundary-value  problem.  The  archetypal  free- 
space  diffraction  problem  Is  to  determine  the  field  over  the 

I 1 

\ plane  z»L,  for  a wave  propagating  nominally  In  the  +z  direction, 

• given  the  field  over  the  plane  z»0.  Let  (x, z,t ) , where 

? “ (x,y,z),  X ■ (x,y),  be  the  positive  frequency  electric  field 


In  other  words,  the  wave  under  consideration  has  no  components  In 
Its  angular  spectrum  which  transport  power  In  the  -z  direction. 
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operator  for  a scalar  wave  that  propagates  nominally  In  the  +z 
direction  In  vacuum.  The  corresponding  quantum  diffraction  problem 
is  as  follows.  Given  the  quantum  characteristics  of 

2 E^'*’^(x,0,t)  , (2.1) 

determine  the  quantum  characteristics  of 

= E^^^(x,L,t)  , (2.2) 

for  an  arbitrary  positive  L.  Herein  we  show  how  to  formulate 
these  planar  characterizations;  the  propagation  question  is 
addressed  In  Section  II. 2. 

The  classical  positive-frequency  electric  field,  ^^^^(r,t), 
can  be  expanded  on  a fixed  transverse  plane 

^ ■ {(x,y,z):-“<x,y<'»}  , 

and  an  arbitrary  time  Interval  T in  terms  of  a complete  orthonormal 

z 

(CON)  set  of  spatio-temporal  mode  functions  {5^(x,t):  (x,t ) ei^xT^ } , l.e. 

E ^"^^x.t)  H E^'^^J.z.t)  - le  £ (x,t)  , (2.3) 

.z  - « 


where 


As  in  classical  optics,  the  scalar  wave  analysis  Is  not  completely 
rigorous.  However,  because  the  generalization  of  our  quantum 
propagation  theory  to  vector  waves  involves  only  added  notatlonal  burden 
we  restrict  ourselves  to  the  scalar  case  • 


,t)C*(5,t)d5dt-«^,, 


z z 


(2.U) 


n 


and  {e^}  are  the  classical  mode-amplitudes.  When  this  classical 

' field  Is  quantized,  the  classical  mode-amplitudes  are  changed  to 

a set  of  quantum  operators,  {A  a } , where  {A„}  Is  a set  of  mode- 
• n n ^ 

dependent  c-numbers,  and  is  s set  of  modal  photon  annihila- 

tion operators  obeying  the  following  canonical  commutation  rules 


(OCR) 


(2.5) 


The  proper  way  to  carry  out  this  quantization  procedure  Is  not 
described  In  this  paper;  Its  relationship  to  the  usual  volume 
quantization  procedure  will  be  discussed  In  a moment. 

In  general,  we  choose  not  to  work  with  the  electric  field 
operator  obtained  by  the  preceding  method,  but  rather  with  the 
field  operator 


n 


(2.6) 


obtained  by  dropping  the  X^.  For  quaslmonochromatlc  (narrowband) 
fields  In  free  space  with  a nominal  radian  frequency  the  field 

(2.6)  Is  proportional  to  the  positive-frequency  electric  field 


operator,  viz 


E^^‘*‘^(x,t)  - (ftw^/2€Q)^'^^E^(x,t)  , 


(2.7) 


i 

i 

I 

I where  Is  the  vacuum  permittivity.  The  quaslmonochromatlc 

I 

I condition  Is  not  actually  required  In  our  theory  once  It  Is  under- 

I stood  that  we  are  working  with  the  field  E , rather  than  the  usual 

electric  field  E For  simplicity,  however,  quasimonochromaticity 

will  be  assumed  henceforth. 

Prom  (2.^)  and  (2.6)  we  have  the  field  commutator 

[E  (x,t)  , E^(x’,tM]*  6(x-x’ )fi(t-t ' ) , (2.8) 

where  E]t(3f,t),  the  adjoint  operator  of  E (x,t).  Is  the  negative- 
Z z 

frequency  counterpart  of  the  field  operator  E (x,t).  Equation 

z 

(2.8)  Imposes  a condition  to  be  fulfilled  by  the  quantized  field  at 
any  transverse  plane,  and  forms  the  foundation  of  our  propagation 
analysis. Apart  from  a proportionality  constant,  Helstrom  has  shown 
[4,  Ch,5]  that  (2.8)  gives  the  electrlc-fleld  commutator  in  the 
limit  of  paraxial  propagation.  More  generally,  by  excluding 
evanescent  waves  and  augmenting  (2.6)  to  Include  waves  which  pro- 
pagate nominally  In  the  -z  direction,  we  can  obtain  the  usual 
volume-quantization  commutator  [4],  [10],  [15]  for  the  electrlc- 
fleld  operator  E from  a planar-quantization  procedure,  and 

^Z 

vice  versa.  We  shall,  however,  purposefully  retain  the  evanescent- 
j wave  contribution  to  (2.6)  because  of  Its  role  In  preserving  the 

field  commutator  (2.8)  at  all  planes  transverse  to  the  nominal 
direction  of  wave  propagation.  Indeed,  the  validity  of  (2.8) 
at  any  z-plane  will  force  us  to  Identify  the  quantum-noise  sources 
responsible  for  Its  preservation,  thus  making  evident  to  what 
extent  we  can  control  the  quantum  state  at  the  receiver  In  subsequent 
■'  communication  analyses  (cf.  Sections  IV,  V). 
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A complete  specification  of  the  quantized  electromapinetlc 
field  requires.  In  addition  to  (2.i*),  (2.6),  (2.8),  knowledge 
of  the  density  operator,  p,  describing  the  states  of  the  modes. 

In  the  absence  of  classical  statistical  dependence,  we  have  that 

P ■ « Cj,  , (2.9) 

n 

where  p describes  the  state  of  the  n^^  mode.  The  state  of  a 
n 

quantized  aperture  field  E.  (x,t),  the  field  E (x,t)  restricted  to 

an  aperture  A In  L , can  be  given  similarly  by  a modal  ex- 
z z 

pension  with  a density  operator  specifying  the  state  of  the  modal 

photon  operators.  If  the  planar  field  E„(x,t)  Is  specified 

z 

first,  the  state  of  E.  (x,t)  can  be  determined  from  that  of  E (x,t) 

*z  * 

by  modal  transformations  (see  Appendix). 

To  solve  the  propagation  problem  we  must  show  how  the 

6-functlon  commutator,  (2.8), Is  preserved  and  how  the  density  operator 

Is  transformed  by  free-space  diffraction. 

II. 2 Modal  Theory  of  Propagation 

Consider  the  llne-of-slght  free  space  optical  communication 
channel,  whose  propagation  geometry  Is  shown  In  Fig.  2.  A 
quasimonochromatlc  quantized  radiation  source  produces  a field 
In  the  z*0  plane  that  Is  propagating  nominally  In  the  +z  direction 
and  la  space-limited  to  a transmitter  exit  aperture  Aq.  After 
propagation  through  L meters  of  free  space,  the  field  is  collected 
over  a receiver  entrance  aperture  A^^.  For  convenience,  we  shall 
assume  that  A^  and  Aj^  are  circular  apertures  of  diameters  d^  and 
dj^  respectively  whose  centers  lie  on  the  z-axls.  The  quantum 
propagation  problem  Is  simply  stated:  given  the  quantum  field 
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j^roduced  by  the  source  over  the  transmitter  aperture,  determine 
the  quantum  field  that  results  over  the  receiver  aperture.  A 
more  precise  formulation  follows. 

Suppose  the  source  produces  a field  In  the  transmitter  pupil 
that  Is  time  limited  to  an  Interval 


t -T<t<t  } . 
U O — “ o 


We  denote  the  field  operator  In  the  transmitter  aperture  by 

E,  (x,t).  This  operator  has  the  mode  expansion 
*0 


(x,t)*  * 


n 


(2.10) 


where  Is  a CON  set  on  AqxTq  , and  the  {a^}  satisfy  the  CCR. 

To  complete  the  quantum  aperture  field  specification,  we  assume 
the  {a_}  In  (2.10)  have  a known  density  operator  p.  . 

n Hg 

In  a similar  manner,  we  assume  that  the  receiver  time  limits 
the  field  In  its  entrance  aperture  to  an  Interval,  Tj^,  that  Is 
sufficient  to  encompass  the  entire  source  transmission.  We 
denote  the  field  operator  In  the  receiver  aperture  by  E.  (x,t), 
which  we  expand  as  follows 


E.  (x.t)-  Ib  c ^(x,t)  . 
*L  n 


(2.11) 


In  Eq.  (2.11)  the  form  a CON  set  on  and,  according  to 
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our  planar  quantization  procedure,  the  must  satisfy  the 

CCR.  To  solve  the  propagation  problem  we  must  determine  p.  , 

the  density  operator  for  the  ^b  ),  from  p.  . Our  solution 

n Aq 

to  this  quantum  field  problem  must  subsume  the  classical  field 
behavior,  i.e.,  the  average  received  field 


(x,t)STr|p^  (x,t))  (2.12) 

li  L L 


must  be  given  In  terms  of  the  average  transmitted  field 


E.  (x,t)STr(p.  E.  (x,t))  (2.13) 

*0  0 0 


by  the  classical  diffraction  formula  . 

Conceptually,  it  Is  convenient  to  solve  this  propagation 
problem  by  a three-step  procedure: 

1)  determine  from  E.  (x,t)  and  P.  the  quantum  characteristics 
of  EQ(x,t),  the  field  operator  over  the  entire  z»0  plane, 
l.e.,  for  (x,t)  c 

11)  determine  from  EQ(x,t)  and  Its  density  operator,  Pq, 
the  quantum  characteristics  of  Ej^(x,t),  the  field  operator 
over  the  entire  z»L  plane,  i.e.,  for  (x»t)  e 
ill)  determine  from  E, (x,t)  and  its  density  operator,  p. , 

ii  M 

the  quantum  characteristics  of  E.  (x,t). 

Within  the  foregoing  schema,  tasks  1)  and  111)  are  trivial. 
Specifically,  we  can  decompose  E^(x,t)  Into  E,  (x,t)  plus  a 

U Aq 

field  operator  defined  on  (Eq-Aq)  * Tq.  All  the  spatio-temporal  modes 
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of  the  latter  operator  must  be  in  the  vacuum  state  because  we  assume 
that  the  source  acts  on  the  half  space  z^O  only  through  the  aperture 
Aq,  and  we  assume  that  there  is  no  classical  interference  present, 
such  as  thermal  noise.  Thus  if  we  exnand  EQ(x,t)  ii  the  form 

EQ(x,t)*  * (2.14) 

n 

where  {c°}  are  CON  on  T , and  the  {c  } satisfy  the  CCR, 

n u u u n 

we  find  that  task  1)  is  completed  teiklng 

Pn"P.  ®|0><0|  . (2.15) 

^ ^0  n 

Task  ill)  is  conceptually  even  simpler;  it  corresponds  to  ex- 
tracting from  a quantum  planar  field  on  spatio-temporal 

modes  that  span  A^^xT^.  Task  il)  will  be  accomplished  by  exploit- 
ing the  classical  mode  theory  of  diffraction.  The  following 
classical  results  are  well  known  Cl6],Cl7];  they  are  presented  to 
establish  notation  . 

A.  Classical  Fields 

As  in  (2.1),  (2.2  ) let  us  use  EQ(x,t)  and  E^(x,t) 

to  denote  the  classical  fields  associated  with  £Q(x,t)  and  E^(x,t) 
respectively.  In  the  quaslmonochromatic  limit  in  which  we  are 
working  Eq  and  E^  can  be  uniquely  decomposed  into  their  propa- 
gating and  evanescent  wave  components,  viz 

EQ(x,t)-EP(x,t)+EQ(x,t)  , (2.16) 

and  EL(x,t)-EP(x,t)+E  ®(x,t)  , (2.17) 
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where  the  two-dimensional  spatial  Fourier  transforms  of  the 
propagatlng-wave  (e^)  components  vanish  for  spatial  frequencies 
beyond  uJq/2tiC  In  magnitude'  and  those  of  the  evanescent-wave  (E®) 
components  vanish  for  spatial  frequencies  below  u)^/2irc  in  magnitude. 

The  propagating- wave  component  of  E^^  Is  generated  by  the 
propagatlng-wave  component  of  Eq  through  the  quaslmonochromatlc 
diffraction  Integral 


:P(x,t)-  df  f dx«  EP(x',t)GL(X,t;X’,t’ 


) , 


(2.18) 


0 “0 


where  the  Green's  functloh  Xi^' satisfies 


aL(x,t;x' ,t' )-  I df  clrc  ( ) 6(t-t '-L/c(l- | 


‘exp{  -J2TrCLx  '^(1-  1 X.T  l2)^''2-r-(x-X’ ) ] } 


(2.19) 


In  (2.19),  Xq  Is  the  wavelength  associated  with  and  f Is  a 
two-dimensional  spatial  frequency.  The  evanescent-wave  component 
of  will  be  nearly  zero,  regardless  of  E^,  for  the  usual 
case  of  L>>  X^ 

Inasmuch  as  E^  Is  the  restriction  of  E^^  to  and  the 

L 

boundary  condition  on  z-0  is  EQ(x,t)H0  for  (x,t);C AqxTq,  the  foregoing 
results  Imply  that 


E.  (x,t)- 


fdt'  fdx*  E.  (x’,t’)G.(x,t,;r',t')  . 

/o  ^ 


(2.^0) 


Now,  because  is  a compact  operator  [17],  [18]  that  maps 
2 2 

L (AqXTq)  Into  L via  (2,20)  , It  possesses  a modal 

decomposition  of  the  form  ^ 


. (2.21) 

n 

where  and  are  CON  sets,  so-called  output  and  input 

n n 

eigenfunctions,  on  and  Aq’^Tq  respectively.  These  eigen- 

functions are  determined  from  the  Fredholm  equations 


L *^^2  1^2  ^L^^l*^l**2*^2^*’n^*2*^2^“^n^n^^l*^l^  * (222) 

' 0 *0 

f dt*  [dx'i:°(x‘,t')G.(x,tjx',t')-ny^cWx,t)  , (223) 

Ja  n L n n 


where  the  kernel  Is 


^L^*l*^l**2’^2^' 


Jl  (x’,t';X^,ti)GL(x' ,t’ ;X2,t2)  . (2.24) 


^ Strictly  speaking,  to  make  G^^  a compact  operator  we  should 
replace  the  6-function  In  (2.19)  l^y  the  bandlimited  6-functlon 
W slnc[W(t-t'  - L/c(l-| X^f I ^)^^^)]  for  some  finite  W.  The  parameter 
W is  the  bandwidth  of  the  quaslmonochromatlc  radiation.  As  the  proof 
of  the  quantum  propagation  theorem  and  all  subsequent  manipulations 
can  be  carried  out  with  some  additional  labor  using  the  bandlimited  Gj^, 
we  have  chosen  to  present  the  simpler  formulation. 
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The  eigenvalue,  n^,  associated  with  nth  Input-output  eigen- 
function pair  is,  physically,  the  fractional  energy  transfer  from 
Aq  to  results  when  we  transmit  ?j^^(x,t).  Because 

the  free-space  channel  is  a passive  non-absorbing  propagation 
medium,  we  have 

0<  n^<l  , for  all  n.  (2.25) 

There  are  two  loss  mechanisms  that  make  n„<l.  First,  there  is 
the  radiative  loss  that  is  incurred  when  Aj^  is  not  large  enough 
to  encompass  the  entire  z*L  plane  field  that  results  from  trans- 
mission of  5^^(x,t).  Second,  there  is  the  evanescence  loss  that 
n 

is  incurred  because  5yj^(x,t)  contains  spatial-frequency  components 
beyond  in  magnitude. 

By  use  of  the  mode  decomposition  for  0^  we  can  replace  the 
diffraction  integral  (2.20)  with  an  equivalent  parallel-channel 
(modal)  description.  Specifically,  if  we  expand  e . using  the 

input  eigenfunctions  as  the  basis  set,  viz 


(x,t)-  (2.26) 
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we  find  that  E.  has  the  expansion 

Ea  (i.t)-  (2.27) 

L n 

In  terms  of  the  output  eigenfunctions  • 

B.  Quantum  Fields 

Just  as  In  the  classical  case,  let  us  decompose  the 
field  operators  Eq  euid  Into  their  propagating  and  evanescent 
wave  components,  l.e. 

EQ(x,t)  - EQP(5,t)  + EQ®(x,t)  , (2.28) 

EL(x,t)  - ELP(x,t)  + EL®(x,t)  . (2.29) 

We  are  Interested  In  developing  the  quantum  characteristics  of 
E^  given  those  of  Eq.  From  the  classical  description  of  the 
propagation  process.  It  Is  apparent  that  E®  Is  not  Influenced 

ii 

by  En  In  the  usual  case  of  L»1  . Thus,  the  annihilation  op- 
u o 

erators  In  any  mode  expansion  of  E^  must  all  be  In  the  vacuum 

P 


state,  and  we  can  concentrate  our  attention  on  E, 


(2.32) 


L n 


In  terms  of  the  classical  ouput  eigenfunctions,  then  the  photon 
operators  {b^}  satisfy 


(2.33) 


Where  {n^j)  are  the  classical  eigenvalues  amd  {c^}  are  vacuvun- 
state  photon  operators.  The  density  operator  for  the  {b^> 
Is  determined  from  p,  and  (2.33)  through  iA.2),  (A.l|).  In 
particular.  If 


P.  ■ ®p. 

^0  n ^Qn 


(2.3^) 


we  obtain 


V f V 


(2.35) 


where  for  each  n p^  Is  computed  from  p^  ^ via  (A.9). 

L 0 


Proof: 


We  have,  from  (2. 28)- (2. 30)  and  the  surrounding  discussion. 


that 


Ej^(x,t)« 


|dt’  jdx'EQ(x',t')GL(x,t;x't»)+EL®(x,t) 


. (2.36) 


, ! 

t \ 


t 


Restricting  (2.36)  to  and  using  the  fact  that  E^Cx.t)- 


(x,t)  on  Aq’^^q  we  find 


E^  (ic,t)-  |dt’  I 


dt’  |dx'E^^(x;t)Oj^(x,t;x' ,t* ) 


^ ^0 


(2.37) 


+ J Jdt 'dx'EQ(x’ ,t' )Gj^(x,t;x' ,t' ) + Ej^®(x,t)  . 

(ZqXT)-(AoXTo) 

We  now  use  (2.37)  to  evaluate  the  series  (2.32)  and  obtain  with 
the  aid  of  (2.21) 


where  d^  Is  the  operator 


(2.38) 


f I 

Jdt  |dxi:^^*(x.t)  ! ||dt’dx'EQ(x;t)GL(x,t;x‘t»)+EL®(x,t) 


(ZqXT).(AoXTo) 


•(2.3i 


Because  {b„)  obey  the  CCR,  It  follows  that 
n 


(2.il0) 
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In  terms  of  a set  of  photon  operators  {c^}  and  arbitrary  phase 

factors  Furthermore,  since  EQ(x',t')  la  a vacuum-state 

field  on  (ZqxT )-(AqXTq)  as  Is  Ej^®(x,t)  on  A^xT^,  the  are  all 

vacuum-state  operators.  For  any  {41^}  the  state  of  the  modes 

described  by  the  photon  operators  {exp(J^^)c^}  is  therefore  the 

vacuum  €)|0><0|  , so  (2.33)  and  the  theorem  follow  immediately  with 
n 

the  Inconsequential  choice  the  density  operator  transformation 

is  given  In  the  Appendix. 

Remarks : 

The  basic  physical  content  of  the  quantum  propagation  theorem 

is  as  follows.  For  and  E.  expanded  In  the  series  (2.31)  , 

*0  *L 

(2.32)  with 


“n  • 


(2.41) 


the  classical  mode  theory  of  diffraction  Implies  that 


(2.42) 


Now,  If  n^<l  we  cannot  have 
n 


b -an  , 

n n n 


because  b„  would  violate  the  OCR.  There  has  to  be.  In  this  case, 
a photon  operator  c^,  that  satisfies  the  CCR,  which  contributes 
to  bj^  according  to  (2.33).  Physically,  when  ^ 1 the  output 

eigenfunction  5 ^^(x,t)  couples  to  modes  that  are  not  generated  by 
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E.  (x,t).  Because  all  such  fields  are  In  vacuum  state,  (2,33) 
is  consistent  with  the  classical  mode  result  (2.42).  The  develop- 
ment In  Section  V provides  a detailed  example  of  the, role. such 


auxiliary  modes  play  In  preserving  the  commutator  (2.8),  as  well 
as  the  effect  of  auxiliary  modes  on  the  extent  to  which  the  receiver 
state  can  be  controlled. 

The  quantum  propagation  theory  developed  In  this  section  is 
actually  more  general  than  the  preceding  theorem  Indicates.  In 
particular.  It  can  be  applied  to  arbitrary  linear  propagation 
media,  and  can  handle  the  case  of  non-vacuum  c^  In  (2,33) » 
c.f.  Section  V.l. 

C,  Coherent  State  Propagation 

An  ordinary  optical  source  Is  one  whose  radiation  field  Is 
in  a Glauber  coherent  state  or  a classically  random  superposition 
of  such  states.  In  free-spece  quantum  receiver  theory,  quantiza- 
tion of  such  a source  Is  generally  performed  (conceptually)  In 
the  receiver  entrance  aperture  or  fictitious  measurement  cavity, 
even  though  the  field  Is  physically  quantum  mechanical  In  the 
neighborhood  of  the  source.  By  means  of  the  quantum  propagation 
theorem  we  can  rigorously  Justify  the  use  of  classical  propagation 

theory  for  ordinary  sources. 

Suppose  that  the  quantized  transmitter  field,  E.  (x,t), 

has  a density  operator  In  the  modal  expansion  (2.13)  whose  P- 
representatlon  possesses  a classical  analog,  l.e..  It  has  the 
coherent  state  expansion 


)• 


o ><a 
n n 


)P(aj^,a2, 


n 


)nd^aj^A 


(2.4 
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where 


Is  a classical  probability  density.  The 


P niA 
n 

quantum  propagation  theorem  (2.31)-(2.33)  In  conjunction  with  the 
density-operator  transformation  calculus  presented  In  the  Appendix 


can  be  readily  used  to  show  that  the  receiver  field,  E.  (x,t), 
has  the  density  operator 


((3|n  I )P(a 

- n n 'n  n' ' ' 
n 


1*2 


• • • 


n 


)nd2ayir  ^2, 


In  the  modal  expansion  (2.32)^.  Equations  (2.^3>  (2.^^)  demonstrate 
that  the  received  field  behaves,  quantum  mechanically,  as  though 
the  propagation  were  classical  eind  quantization  were  applied  at  the 
receiver.  We  shall  see.  In  the  sequel,  that  the  quantum  propagation 
theorem  has  slgnlglcant  content  for  more  general  quantum  transmitters. 


Equation  (2.44)  has 
Interference  or  noise. 


to  be  modified  In  the  presence  of  classical 
which  Is  not  our  concern  here. 


Ill  Two-Photon  Coherent  States 


We  now  Introduce  a class  of  radiation  states  that  are  closely 
related  to  ordinary  coherent  states,  but  which  possess  a number 
of  advantages  from  a communication  viewpoint.  A comprehensive 
discussion  of  the  mathematical  and  physical  nature  of  these  states 
has  been  given  elsewhere[7] « so  that  the  present  discussion  serves 
primarily  to  establish  notation  and  Indicate  the  possible  comm- 
unication Implications. 

III.l  State  Specification  and  Generation 

For  a single  radiation  mode  of  frequency  with  photon 
annihilation  operator  a,  the  two-photon  coherent  states (TCS) 

s 

lS>g  are  defined  to  be  the  eigenstates  of  a transformed  photon-  • - - 

operator  b with  eigenvalues  B.  Specifically,  we  define 

' 

t ' 

b«ua+va,  (3.1) 

for  complex u andv  such  that 


the  TCS  I satisfies 
S 


(3.2) 


b|B>g-  6|B>g  . (3.3)  j 

( 

When  v«0,|6>  becomes  the  ordinary  coherent  state  |S>. 

O 

For  non-zero  v,  the  state  p»|6>  <6|  does  not  possess  a well- 

oo 

behaved  P representation. 
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The  TCS|i>>  has  the  following  important  fluctuation  behavior 

s 


o.<6|a|6>^  - u«e-  v6»50  , 
8 g 


g<B|a'‘’a|e>g  - |6|^  + |v|2  , 


(3.4) 


(3.5) 


<Aaj^^>  - |u-v|Vn  , <L^>  - |u  +v|2/4  , (3.6) 


where 


asa,  + Ja, 


(3.7) 


for  self-adjoint  a^^,  a.2  defines  the  quadrature  components  of  the 
mode  in  question  and  Aaj^=a^-<aj|^>  is  the  fluctuation  in  a^,  1»1,2, 
(In  these  expressions,  the  angular  brackets  denote  quantum  ex- 
pectation with  respect  to  a given  density  operator.) 

For  arbitrary  u,  v,  the  TCS|B>  achieves  the  minimum  uncer- 

O 

talnty  product 


<Aa'^^><Aa^^>  - Z/l6  , 


(3.8) 


for  the  quadrature  components  of  the  photon  operator 


a'=a  exp  J J tan 


-1 


.1(u  V - V u) 


2|ii  1 1 v|  + w*v  ♦ v*u 


(3.9) 
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with 

<Aai^>  - (|w|  - |v|)^/i»  , <Aa^2>-  (|m|  + |v|)^/4  .(3.10) 

When  p/v  is  real,  TCS  are  the  usual  minimum  uncertainty  states  | 

I for  the  Heisenberg  lnequallty<Aaj^‘^><Aa2‘^>il/l6.  In  the  limit  | 

i 

W»v**«  with  u/v-^1  and  u/\h— 1,  TCS  become  the  eigenstates  of  a^^  and  | 

a2  respectively,  viz  . 

4 

2 

I <Aa^  >-^0  , for  p/v-^1  , 

(3.11) 

2 ^ 
<Aa2  >'►0  , for  p/v-^-l  . 

It  follows  from  (3.6)  or  (3.10)  that  the  quantum  fluctuations 
In  a^  and  a2  (or  a^  and  ap  can  be  exchanged,  at  constant  min- 
imum uncez^talnty  product,  by  V2u*ying  y,v  In  a TCS.  For  coherent- 

2 2 

state  radiation  we  have  that  <Aa^  > « <Aa2  >■  1/4  , and  random  super- 
position of  coherent  states  always  Increases  the  fluctuation 
' levels.  Thus,  the  exchange  of  quadrature-fluctuation  levels  Is 
not  a feature  of  conventional  sources. 

Whereas  a coherent -state  radiation  field  Is  generated  by  a 
one-photon  laser  operated  well  above  threshold,  a TCS  can  be 
I generated.  In  principle,  by  a two-photon  laser  operated  well 

above,  threshold.  The  theoretical  relationship  between  TCS  radla- 

^ • 

tlon  and  stimulated  two-photon  processes  has  been  elaborated  In  [7]. 

As  yet,  no  TCS  radiation  has  been  experimentally  observed,  although 
the  generation  of  such  light  remains  a physical  possibility. 
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Naively,  the  rationale  for  improving  free-space  optical 


1 


communication  performance  through  use  of  TCS  radiation  is  as  follows. 

At  the  transmitter,  we  impress  analog  or  digital  information 
onto  the  average  value  of  the  a^^  quadrature  of  a single  TCS  radia- 
tion mode.  By  controlling  u,v  we  can  reduce  the  quantum  fluc- 
tuation on  the  signal  (a^)  quadrature  below  the  level  for  co- 
herent states.  Thus,  measurement  of  at  the  receiver  yields  a 
higher  slgnal-to-nolse  ratio  for  analog  communication  or  lower 
error  probability  for  digital  communication  than  that  obtained 
with  a coherent-state  transmitter  of  the  same  energy. 

1 

i 

III. 2 State  Propagation  j 

I 

The  preceding  description  of  potential  TCS  communication  1 

t 

advantages  is  oversimplified  because  it  neglects  the  effects  of 

propagation  from  the  transmitter  to  the  receiver.  Consider  the 

line-of-sight  free  space  propagation  channel  that  was  analyzed  1 

in  Section  II,  Suppose  that  a TCS  transmitter  produces  a field  \ 

(x,t)  whose  state,  , in  the  expansion  (2,31)  is  j 

0 *0  1 


Pa  “ |8>m<3|®|0><0| 
“0  °°  n 


(3.12) 


i.e.,  the  first  input  eigenfunction  is  in  a TCS  and  all  the  re- 
maining eigenfunctions  are  in  vacuum  state.  Via  the  quantum  propaga- 
tion theorem  we  can  separate  the  state  generation  and  propagation 
problems  and  find  that  E.  (x,t),  the  receiver  field,  has  state 

Pa  -Pj^®|0><0|  (3.13) 

L n 


i 
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In  the  expansion  (2.32).  in  (3.13) > Is  the  state  of  the  first 
output  eigenfunction;  It  Is  the  state  of  a photon  operator 

b - a + (3.14) 

where  a,  c are  photon  operators  that  are  In  states  |&>  an-l  1 0> 

O 

respectively.  The  output  field  state  can  be  found  from  (A. 11), 
and  Is  In  general  a mixed  state. 

For  the  purpose  of  eliciting  the  effect  of  propagation  on 
TCS  quantiim  noise  reduction.  It  Is  simpler  to  calculate 

the  fluctuation  behavior  of  b directly  from  (3.4)- (3. 7)  and  (3*14)  than 
to  use  the  transformed  state  (A. 11).  We  find  that 


<b>  « , 

(3.15) 

<b‘‘'b>  - n^dBl^  + lv| 

1^)  . 

(3.16) 

<Abj^^>  ■ + 

(1-t1i)]/4  , 

(3.17) 

<Ab2^>  ■ + vj^ 

+ (i-ni)]/4  , 

(3.18) 

where  b=b^  + Jb2  defines  the  quadrature  components  of  b. 

Comparison  of  (3-6)  and  (3.18)  reveals  the  following  behavior: 

p 

the  TCS  noise  reduction  that  Is  manifest  on  a^^  when|y-v|  < 1 Is 
nearly  absent  from  b^  unless  Physically,  when  n2^<<l 


the  quantum  characteristics  of  b are  dominated  by  the  vacuum  co- 
herent state  properties  of  c through  (3.1^) • Thus,  to  reap  the 
communication  benefits  suggested  by  (3.6)  we  should  arrange  the 
propagation  geometry  so  as  to  achieve  a near  unity  fractional 
energy  transfer  (eigenvalue)  for  the  TCS  transmitter  mode;  such  a 
situation  will  prevail  In  the  classical  near-field  propagation 
regime  [19].  Alternatively,  for  the  classical  far-fleld  propagation 
regime  in  whlchn^<<l,  we  could  conceivably  Improve  on  (3.17)  by 
use  of  a TCS  local  oscillator  near  the  receiver  to  cause  the 
photon  operator  c In  (3.1**)  to  be  in  state  1 0>  . These  near  field 

O 

and  far  field  cases  are  considered  In  Sections  IV  and  V respectively. 


IV  Quantum  Transmitter  Optimization 


( 

t 

i 

i 

I 

I ' 

I 


The  discussion  surrounding  (3*l^)-(3*l8)  has  the  following 

generalization.  Suppose  a quantum  transmitter  produces  a field 

operator  E.  (x,t)  with  state 
*0 


Pa  ■ |><|®l0><0|  , (il.l) 

*0  n 

where  )>  Is  an  arbitrary  pure  state.  When  nj^<<l  the  resulting 
state,  p^,  for  the  first  output  eigenfunction  In  the  expansion  of 
E.  (x,t)  will  be  approximately  a coherent  state  (Insofar  as  field 
fluctuations  are  concerned) 


Pl^hi^'^^  <a»<ni^^^<a>| 


(i».2) 


In  the  absence  of  receiver  state  control,  where  a Is  the  photon 
operator  associated  with  the  first  transmitter  mode.  Evidently, 
with  a field-measurement  receiver,  whatever  benefits  may  be 
derived  from  optimum  quantum  state  generation  at  the  transmitter, 
they  are  limited  to  the  regime  wherein  n^^-l  when  no  receiver 
state  control  can  be  employed.  It  therefore  behooves  us  to 
briefly  address  the  eigenvalue  behavior  of  the  classical  free- 
space  mode  decomposition  before  attempting  a quantum  transmitter 
theory. 


IV. 1 Near  Field  and  Far  Field  Propagation 


For  the  usual  case  in  which  the  aperture  diameters  dQ  and  dj^, 
and  the  path  length  L,  are  all  much  larger  than  the  nominal  wave- 
length  evanescence  loss  can  be  neglected  and  hence  the  Fresnel 
approximation 


Q^(«,tjx’t')  » 6(t-t'+L/c) 


(-JXoL)'^  exp[-J2itL(l+|x-x' |^/2L^)/X^]  , 


(4.3) 


can  be  used  in  place  of  (2.19)  • The  mode  decomposition  (2.21)  - 
(2.24)  associated  with  the  Fresnel  approximation  Green's  function 
is  well  known  [19]i  [20];  the  salient  features  are  summarized 


below. 


Because  of  the  delta  function  time  dependence  in  ( 4.3  ), 
if  we  teike  » Tq-L/c  we  obtain  the  factored  mode  decomposition 
{C^°(x)Xnj(t) ; 5n^(*)Xj„(t-L/c)  ;n^}  where  the  temporal  input  modes 
{Xjjj(t)}  are  an  arbitrary  CON  set  on  Tq  and  the  spatial  mode  de- 
composition satisfies 

|dX2(dj^J^(Trdj^|  Xg-x^l  /\^L)/2\^L\  X2-x^|  )exp[jTT(  |xj^|  ^-|x2l  ^)/XqL] 


Cn  ^*2^  " ’^n^n  ^^1^  * 


(4.4) 


I dx'j;jj°(x')(-JXQL)"^exp[-J2TrL(l+j  jc-x'  |^/2L^)/X^] 

An 


(11.5) 


F 


'I 


Explicit  results  for  t^j^^Cx);  C^^^Cx);  n^)  are  available  [20]. 

The  eigenvalues  depend  parametrically  on  the  Fresnel  number 

D - ( (M.6) 

and  the  input  and  output  eigenfunctions  are  prolate  spheroidal 
wavef unctions . For  our  purposes  it  is  sufficient  to  note  the 
asymptotic  behavior  of  the  eigenvalue  distribution.  Suppose  the 
eigenvalues  are  arranged  in  decreasing  order,  then  in  the  far- 
fleld  region,  D<<1,  the  maximum  eigenvalue,  is  approximately  equal  to 

D,  and  all  other  eigenvalues  are  insignificant.  When  D>>1,  near-field 

propagation  prevails;  there  are.Jthen  p eigenvaluea  near,  1 and  the 
rest  are  insignificant. 

IV. 2 Near-Field  Transmitter  Optimization 

Consider  the  following  problem  of  quantum  state  optimization 
in  near-field  free  space  communications.  A continuous  real- 
valued random  variable  a with  density  function  p(a  ) is  to  be 
transmitted  from  AqxT^  by  placing  the  photon  annihilation  operator 
a for  the  mode  Cj^^(x)x^(t)  in  state  p . Let  us  assume  that  llneau* 
modulation  is  employed  so  that  when  a is  chosen  for  transmission  , 
is  constrained  to  satisfy 

<a>=Tr(Pjja)«aK  , (4,7) 

where  K is  a positive  constant.  The  receiver  measures  b^^,  the 


-31- 


cosine  quadrature  of  the  photon  annihilation  operator  b associated 
with  the  mode  x)X2(t-L/c)  on  A^xT^.  Because  we  are  in  the 
near-field  propagation  regime,  we  may  assume  n^~l.  and  hence 
b*a,  bj^'a^^.  The  problem  is  to  choose  p^and  K to  maximize  the 
average  signal-to-nolse  ratio 


SNR 


a 


av 


T ? 

dop(o)Tr[Pjj(Aaj^)  ■ 


(4.8) 


where  a is  the  mean-square  value  of  the  parameter  a,  subject  to 
the  average  photon  number  (energy)  constraint 

|dap(a)Tr(p^a^a)<Ng  , (4.9) 

on  the  transmitter  state.  The  classical  average  transmitted 
energy  is 

|dap(a)|Tr(p^a)|^»K^^  . (4.10) 

Thus,  because<a^a>>j<a>l  (4.9)  is  a more  severe  and  more  appro- 
priate energy  constraint  for  the  quantum  transmitter  problem  than 

2~^ 

the  classical  constraint  K a <Ng. 

We  now  show  that  the  optimum  state  and  modulation  con- 
stant K°  satisfy: 


p °-|k°o(2N„+1)^^^»  <K°o(2N  +1)^''^!  , (4.11) 

a S gg  a 

K°-[o2(2N,+l)/NWg +1)]“^'^^  , 

B 8 ^ 


(4.12) 


where  the  TCS  1 K®a{2N,+l)^''^>  has 

s g 


p-(Ng  + l)/(2Ng+l)^/2  ^ V-Ng/<2Ng  + l)^''^ 


(4.13) 


That  the  optimization  depends  on  p(a)  only  through  o should  not 


be  surprising  because  SNR„  is  a second-moment  quantity.  The 

a V 


proof  that  (4. 11)- (4. 13)  constitute  optimum  choices  of  and  K 


Is  as  follows. 


By  means  of  (3.4)-(3.6)  It  Is  easy  to  show  that  (4.11)-(4.13) 


satisfy  (4.7),  (4.9)  with  equality,  and  give 


SNR^/.HN^CN^+l)  , 


(H.IH) 


Which  we  must  show  to  be  the  optimum  SWR  value.  We  first  re- 

S V 


write  the  energy  constraint  (4.9)  In  the  equivalent  form 


Tr[p(Aa^)^]+Tr[p(Aa2)^]+K^^Ng+l/2  , 


(4.15) 


where 


p=  |dop(a)p^  . 


(4.16) 


A simple  convexity  argument  shows  that  p ° and  K°  will  satisfy 


(4.15)  with  equality.  Moreover,  when  (4.15)  Is  an  equality 


(4.8)  can  be  expressed  as 


SNRg^-  (Ng+l/2-Tr  [p  ( ASj ) ^ ] )/Tr  [p  ( Aa^^ ) ^ ]-l 


(4.17) 
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Now,  because  p is  a density  operator  we  have  the  uncertainty  rela 
tlon 

TrLP(Aa^)^]*Tr[p(Aa2)^]il/l6  , (^.l8) 

and,  from  (4.17),  It  Is  clear  thct  will  satisfy  (4.l8)  with 
equality.  Under  this  minimum  uncertainty  condition  we  have 

SNR  -(N  +l/2)/Tr[p(Aa, )^]-(4TrCp(Aa,)^])"^-l  , (4.19) 
av  s J.  J. 

which  has  the  unconstrained  maximum 

When 

Tr['p(Aa^)^-C4(2Ng+l)]‘'^  . (4.21) 

The  proof  Is  now  complete  In  that  (4. 11)- (4. 13)  satisfy  the 
constraints  and  achieve  (4.19).  Note  that  the  use  of  ordinary 
coherent-state  radiation  In  this  linear  modulation  problem 
results  In  the  maximized  slgnal-to-nolse  ratio 

SNR^^-4Ng  , (4.22) 

when  K-(Ng/o^)^^^  , Pj-lKaXRaj  . 

For  N>>1,  the  optimum  TCS  performance  Is  vastly  superior  to  the 
s 

conventional-source  result. 


(4.20)  ! 
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Note  that  we  have  not  performed  a combined  transmitter- 


receiver  optimization;  the  measurement  of  a^^  Is  fixed  In  advance. 
For  linear  modulation  of  TCS  fields,  however,  one  can  easily  show 
that  the  a^^  measurement  actually  minimizes  the  measurement 
variance  for  any  Input  signal  parameter.  Thus,  the  performance 
(iJ.20)  cannot  be  further  Improved  by  receiver  optimization. 

The  possibility  Is  open  that  a Joint  transmitter-receiver 
optimization  will  lead  to  a superior  performance  than  (4.20), 
which  seems  unlikely  In  the  case  of  a continuous  parameter  a, 

i 

f , but  which  Is  obviously  true  In  the  discrete  case  by  transmitting 

photon  eigenstates. 


V Far-Fleld  State  Optimization 


Recall  from  Section  IV  that  under  far-field  propagation 
conditions  the  quantum  field  fluctuation  characteristics  of  (x,t) 
will  be  dominated  by  the  vacuum  coherent-state  behavior 
Introduced  by  the  radiative  loss  in  the  channel.  In  this 
section,  we  shall  indicate  how  a TCS  source  might  be  used  as 
a receiver  local  oscillator  to  achieve  a state  optimization 
analogous  to  that  in  Section  IV, 2.  The  present  discussion 
is  doubly  speculative  in  that  it  assumes  a source  that  is 
as  yet  unrealized  (two-photon  laser)  and  a spatial-mode 
mixing  for  which  a practical  configuration  has  yet  to  be 
found.  However,  the  predicted  performance  increment  is  so 
large  that  continued  implementation  research  is  clearly 
Justified.  In  particular,  we  shall  find  that  by  means  of 
a TCS  local  oscillator  we  cam  increase  post-measurement 
signal-to-nolse  ratio  by  a factor  equal  to  the  radiative 
loss  of  the  channel. 

Suppose,  as  in  Section  IV. 2,  that  a continuous  real- 
valued random  variable  a is  to  be  transmitted  using  linear 
modulation  from  by  placing  the  annihilation  operator 

a for  the  mode  ?^°(x)x2(t)  in  state  p^.  The  receiver 
measures  b^,  the  cosine  quadrature  of  the  annihilation 
operator  b associated  with  the  mode  Cj^^(x)x2^( t-L/ c)  on 
We  shall  assume  far-field  propagation  prevails  so  that 

b-ni^^^a+(l-ni)^^^c  (5.1) 


A 
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where  r\ (w dQd^/4x  ^L)  <<  1,  and  c Is  a photon  annihilation 

operator  whose  state  Is  unaffected  by  . Note  that  any 

1/2 

phase  factor  that  would  have  multiplied  (l-n^^)  c In  (2.1<0) 
has  been  absorbed  Into  the  photon  operator  definition  for  c 
In  (5.1).  Our  objective  Is  to  maximize  the  slgnal-to-nolse  ratio 


c3a  p(a)Tr(pj^b^)‘^ 
da  p(a  )Tr[p  j^(Ab^)  ] 


(5.2) 


subject  to  the  transmitter  photon  number  constraint  (^.9)  by 
choice  of  p^  and  K.  In  (5.2),  p^^  Is  the  state  of  b,  which. 

In  the  absence  of  any  radiation  source  other  than  the  trans- 
mitter that  produced  E.  will  be  approximately  the  coherent  state, 
(c.f.  (i».2)) 


Pi*  |ni^''^Ko><nj_^^^Ka| 


(5.3) 


Under  (5.3)  It  is  trivial  to  compute  that  the  maximum  slgnal- 
to-nolse  ratio 


SNRav*''''i''s 


results  when  we  use  the  coherent-state  p^*lKa><Kal  with  K»(Ng/a  ) ' 

Comparison  of  (5.^)  and  (l|.22)  reveals  the  well  known  SNR 

av 

attenuation  that  attends  free-space'  propagation.  To  ameliorate 
this  attenuation  by  means  of  state  optimization,  we  must 
control  the  state  of  the  photon  operator  c.  The  details  are  worked 
out  below  , 
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V,1  state  Optimization  by  Spatial-Mode  Mixing 

For  the  assumed  transmission  geometry  we  can  expand  the  z“L 
field  operator  in  the  series > ■ 

Ej^(x,t)«a'>g(x)x3^(t-L/c)+c4>j^(x)X3^(t-L/c)+i:ajj^ejjjCx,t)  , (5.5) 

m 

where  on  Z, 


♦g(x)« 


f.  dx’C  °(x»)(-jX  L)"^exp[-J2TrL(l+|x-x’ 1^/2lV^«]  , (5.6) 


♦r(X). 


(5.7) 


and  {'l>gU)X3^(t-L/c) , <l>R(x)Xj^(t-L/ c) , Cj^(x,t)}  form  a CON  set  on 
Note  that  by  virtue  of  (4.5),  the  restriction  of 
the  "signal"  spatial  mode,  ♦g(x)»  to  is  ^2^^^  ^^^^(x).  By 
a simple  calculation  it  follows  from  (5.7)  and  (2.8)  that 
♦r  (x)  Is  the  other  "received"  spatial  mode  when  we  measure 
^^^^(x)  on  l.e.,  a and  c in  (5.5)  are  the  operators  that 
appear  in  (5.1).  Because  these  operators  are  associated 
with  orthogonal  modes  on  Z x T.  we  may  take  them  to  be  quantum- 

ii  il 

mechanically  Independent,  l.e.,  they  are  members  of  a set  of 
photon  operators  obeying  the  CCR.  Further,  we  suppose  that  when  a is 
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to  be  transmitted  the  mode  described  by  a Is  In  a coherent  state 
(Ng/^)^''^^o><(Ns/^)^'^^o|  (generated  by  the  transmitter), 
and  the  mode  described  by  c Is  In  state  p^,  controlled  by  a source 
(the  local  oscillator)  In  the  vicinity  of  z=L.  Our 
objective  Is  to  maximize  the  slgnal-to-nolse  ratio  (5*2) 
subject  to  the  local-oscillator  energy  constraint 

Tr(p^c''’c)<Np  . (5.8) 

The  preceding  state  optimization  is  straightfoward. 

We  directly  compute  from  (5.1),  (5.2)  that 

. (5.9) 

Rewriting  (5.8)  as 

Tr[p^(Ac^)^]+Tr[P^(Ac2)^]+|Tr(p^c)|^<Nj^+l/2  , (5.10) 

the  reasoning  employed  in  Section  IV. 2 leads  us  to  conclude 

that  the  TCS 


P<,°-l0>gg‘0|  , (5.11) 


yields  the  maximum  slgnal-to-nolse  ratio 


4N. 


SNR 


1+C  ( 1+Nj^)  ( )/rT3^ 


(5.13) 


When  Np>>l,  as  will  surely  be  the  case,  (5.13)  becomes 


SNR, 


o . 


1»N. 


av 


l+(l-n^)/4n3^Nj^ 


(5.14) 


Thus,  when  Nj^>>(l-n^)/4r^  we  find  that  the  post-measurement 
slgnal-to-nolse  ratio  for  the  far-fleld  propagation  channel 
with  a TCS  local-oscillator  receiver  reduces  to  the  near- 
fleld  propagation  coherent-state  transmitter  result  (recall  (4.22)) 


= ““a 


(5.15) 


Apparently  use  of  the  TCS  local  oscillator  permits  the 
post-measurement  quantum  noise  to  be  attenuated  to  the  point 
where  the  transmitter  quantum  fluctuations  dominate  the  slgnal- 
to-nolse  ratio.  One  suspects,  therefore,  that  additional 
benefits  must  accrue  If  we  use  a TCS  transmitter  as  well  as 
a TCS  local  oscillator,  and  Indeed  such  Is  the  case.  We  omit 
the  general  result  for  optimizing  the  slgnal-to-nolse  ratio 
(5.2)  by  choice  of  K,  and  under  the  constraints  (4.9)  and 
(5.8)  In  favor  of  the  following  simpler  special  case.  When 


Nj^  satisfies  Np>>(l-nj^) (Ng+l/2)/2nj^  the^ optimum  slgnal- 

to-nolse  ratio 

^5.16) 

Is  realized  by  using  (^.11)-(^.13)  for  and  K°  and  (5.11),  (5.12) 
for  p °.  Note  that  (5.16)  equals  the  optimized  near- field  result 

w 

(4.14). 

V.2  Implementation  of  Spatial-Mode  Mixing 

Having  found  such  a remarkable  performance  Improvement 
Is  obtained  through  the  beneficial  spatial-mode  mixing  on  of 
a far-fleld  propagation  signal  mode  and  a carefully  chosen  TCS 
receiver  mode,  we  naturally  Inquire  whether  the  results 
are  specious.  Let  us  say  that  there  is  some  finite  circular 
region,  A^,  In  the  z»L  plane  which  Includes  and  Is  such  that 

lA^|0g(x)l^dx:i  , (5.17) 

is  the  region  occupied  by  the  signal  field  In  z=L. 
s 

If  our  receiver  measures  the  <|>o(x)  mode  on  A_  rather  than  Its 

s s 

restriction  C^^(x)  to  Aj^,  we  would  not  be  surprised  to  obtain 

near-fleld  signal-to-nolse  ratios,  because  (5.17)  Is  precisely 

the  condition  that  ?^^(x)  propagates  from  Aq  in  z»0  to  A^ 

In  z*L  without  loss.  The  spatial-mode  mixing  analysis  of  Sec- 

\ 

tion  V.l  therefore  Implies  that  the  effect  of  the  TCS  local 
oscillator  is  to  make  the  Aj^  measurement  of  b^equi valent  to 
the  A.  measurement  of  a, . This  Is  a hollow  equivalence  if  to 

S X 


generate  the  TCS  mode  ♦p(x)  on  we  require  an  apparatus  that 

fills  A . We  need  a technique  to  generate  the  desired  field 
s 

on  A from  a device  that  Is  significantly  smaller  than  A_. 
s s 

As-  yet,  we  have  not  found  such  a technique,  but  neither 

have  we  found  a general  proof  of  Its  Impossibility.  The  one 

result  we  have  obtained  Is  that  the  obvious  answer,  a divergent- 

beam  local  oscillator  located  In  front  of  the  receiver,  cannot 

work.  We  present  a brief  treatment  of  the  divergent  beam 

system,  because  It  illustrates  the  difficulties  encountered 

In  trying  to  realize  the  desired  spatial-mode  mixing  on 

Whether  or  not  any  practical  scheme  exists  for  the  TCS  mode 

generation  remains  an  open  and  important  problem. 

The  geometry  of  the  divergent-beam. receiver  system  is 

shown  in  Pig.  3.  A TCS  source  located  In  the  z-L-^  plane  ('^<L) 

radiates  a diverging  spatial  mode  4>£(x)  from  a circular 

aperture  A . that  is  much  smaller  than  A, . To  achieve  the 
L 'Ll 

desired  mode  mixing  on  Aj^,(J)£(x)  must  give  rise  to  the  field 
♦j^(x)  on  Ag  In  the  z=L  plane.  According  to  the  quantum 
propagation  theorem,  this  propagation  from  A^to  A^  must  be 
(nearly)  lossless  or  else  the  TCS  quantum  characteristics 

will  not  be  preserved.  Thus  we  require  as  a necessary  condition 
the  near-field  propagation  condition 

(^rd^dg/4Ao^)^>l  (5.18) 

where  d^  and  d^  are  tne  diameters  of  A^  and  A^.  For  typical 
parameter  values,  (5.18)  does  not  pose  a serious  difficulty. 
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Unfortunately  (5.18)  is  not  a sufficient  condition.  It  does  not 

suffice,  for  our  present  purposes,  that  there  Is  a mode  ♦^(x)  that 

suffers  hardly  any'  no  radiative  loss 'In-propagation ‘from  to 

There  must  be  a mode  ♦^(x)  on  A^  that  yields  without  radiative  los 

a specific  mode,  determined  from  (5*7)  with  <l>_(x)  modified  to 

s 

Include  the  central  obscuration  of  Aj^  by  A^,  on  A^  (see  Fig.  3). 

Simple  recourse  to  Fresnel  diffraction  theory  demonstrates 

that  the  above  problem  is  over-constrained.  The  far-fleld 

propagation  which  governs  the  signal  mode  4>^(x)  gives  rise  to  a 

s 

2 

quadratic  phase  curvature  exp(-jTr| xH/^^L)  in  the  z*L  plane, 
and  this  same  phase  curvature  persists  in  <|>p(x).  The  divergent 
TCS  beam  has  a quadratic  phase  curvature  exp(-Jir  | x ( in 

p 

the  z=L  plane.  Because  l«h  and  dg  have  already  been 

assumed  (the  latter  condition  is  Implicit  in  (5.18))  the  mode 
generated  by  <l>^(x)  on  the  z«L  plane  will  be  almost 
orthogonal  to  ♦^(x). 
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Transformation  of  Quantum  States 


■ ■■  • 

« 

Appendix 

Let  a^*{aj  : be  a finite  set  of  photon  annihilation 

operators  that  obey  the  CCR.  The  density  operator  describes 
the  state  of  these  modes,  so  that  the  state  , for  the  1th  mode  Is 
obtained  by  tracing  P.  over  the  other  modes.  The  following  Joint 
antl-normally  ordered  characteristic  function  Is  always  well  defined 
[7],  [21] 

Xj(Y,Y»)=Tr[P^  exp(2Y^»aj^)exp(-2Yj^a^^)]  , (A.l) 

where  Y Is  a complex  vector.  Through  normal-ordering  techniques 
we  can  derive  the  following  Integral  representation  for  P [7] 

cl 

^a“  I xf(Y,Y*)exp(-£Yj^a^^)exp(EY»j^a^)d^X/ir^  . (A. 2) 

Let  a**{aj  :1^J£N}  be  another  set  of  photon  operators  which 
are  obtained  from  a by  a unitary  matrix  transformation 

a'-Ua  . (A. 3) 

Because  a obeys  the  CCR  and  U Is  unitary,  a'  will  also  obey  the  CCR. 

of 

The  characteristic  function,  X^<Y,Y*)  for  a'  can  be  computed  from  (A.l) 
and  (A. 3)  with  the  result 

xf ' ( Y . Y» ) -x|(  U^i,  ( U+l) » ) . . ( A . l» ) 
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Thus,  the  state  p , of  a'  Is  determined  by  (A. 2)  and  (A. 4). 

The  above  results  can  be  formally  extended  to  the  case  of  an 
infinite -number  of  modes.  ‘ However,*  Iri  actual  appllcatiohs  there  Is 
usually  some  finite  number  of  modes  that  suffices. 


The  quantum  propagation  theorem  (2.33)  requires  the  first 
half  of  the  following  unitary  transformation 


where  b^,  b^,  a^,  c^  all  obey  the  CCR.  Under  the  condition  i 

i 

l>a,c"®V®  ‘>0  . (*-6) 

n n n n 

where  a={aj},  c={Cj}  , we  have  from  (A. 4) 
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(A. 7) 


n “‘’n  ’ 


where  and  Tor  eachn.p^  can  be  computed  from  part  of  the 


two-dimensional  transformation 


• k “ ! 


r“i/ 

L“'j  . 


f a ] 


(A. 8) 


Where  In  this  case,  (A.^»)  yields 


x^(Y,Y»)=xJ(n^/2Y.n^/2Y»)x;((l-n)^/2Y.d-n)^''V)  . (A.9) 


In  terms  of  the  convenient  density  function  (a,o»)i<o|  p|  a> 


for  coherent  states  la>,  we  have  from  (A.9)  and  (A. 2)  the  convolution 


formula 


p^j;Ja,a«) 


I P^"i6,e«)P^"^((l-n)"^/^(a-n^/^e),(l-n)”^/^(a-n^''^3)»)d^3/Tr(i-n)  . 

J oL  C 


As  an  example  of  these  transformation  rules,  consider  the  case 


In  which  a Is  In  the  TCS  |0**_  with  parameters  W_,  v and  c is  in  the 

S cl  8 


TCS  |0>  with  parameters  M , v , Then  b,  from  (A. 8),  (A.9),  has  the 
2 c c 


characteristic  function 


X^(lf,Y«)-  exp(-Y^AY^/2)  , 


where  Y^  « . . 2. 

“2y, 


(A. 11) 


- (1^-.1£)'- 


A - n* 


-u^v  n-p^v  (1-n) 
a a c c 

1 1 

-J  J 


1/2  . 


, (A.13) 


(A.14) 


Equation  (A. 11)  represents  a TCS  with  added  Gaussian  noise  and  so  Is 

a mixed  state  In  general.  It  can  be  shown  that  (A. 11)  becomes  a pure 

« « 

state  If  and  only  If  ^ Is  real. 

a gL  w w 
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Figure  Captions 

Pig.  1 Generic  block  diagram  for  a single  communication  link. 
Pig.  2 Free  space  propagation  geometry. 

Pig.  3 Divergent-beam  local  oscillator  geometry.  TCS  local 

oscillator  radiates  a diverging  beam  from  the  aperture 
In  z*L-d,  where  Is  much  smaller  than  Aj^  and  t is  much 
smaller  tnan  L.  This  local^oscillator  aperture  prevents 


any  signal  light  from  reaching  the  central  portion  of  A. . 
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